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Abstract — Space-Time block codes (STBC) from Orthogonal 
Designs (OD) and Co-ordinate Interleaved Orthogonal Designs 
(CIOD) have been attracting wider attention due to their 
amenability for fast (single-symbol) ML decoding, and full-rate 
with full-rank over quasi-static fading channels. However, these 
codes are instances of single-symbol decodable codes and it is 
natural to ask, if there exist codes other than STBCs form ODs 
and CIODs that allow single-symbol decoding? 

In this paper, the above question is answered in the affirmative 
by characterizing all linear STBCs, that allow single-symbol ML 
decoding (not necessarily full-diversity) over quasi-static fading 
channels-calling them single-symbol decodable designs (SDD). 
The class SDD includes ODs and CIODs as proper subclasses. 
Further, among the SDD, a class of those that offer full-diversity, 
called Full-rank SDD (FSDD) are characterized and classified. 

We then concentrate on square designs and derive the maximal 
rate for square FSDDs using a constructional proof. It follows that 
(i) except for N = 2, square Complex ODs are not maximal rate 
and (ii) square FSDD exist only for 2 and 4 transmit antennas. For 
non-square designs, generalized co-ordinate-interleaved orthogo- 
nal designs (a superset of CIODs) are presented and analyzed. 

Finally, for rapid-fading channels an equivalent matrix channel 
representation is developed, which allows the results of quasi- 
static fading channels to be applied to rapid-fading channels. 
Using this representation we show that for rapid-fading channels 
the rate of single-symbol decodable STBCs are independent of the 
number of transmit antennas and inversely proportional to the 
block-length of the code. Significantly, the CIOD for two transmit 
antennas is the only STBC that is single-symbol decodable over 
both quasi-static and rapid-fading channels. 

Index Terms — Diversity, Fast ML decoding, MIMO, Orthog- 
onal Designs, Space-time block codes. 

I. Introduction 

SINCE the publication of capacity gains of MIMO systems 
[1], [2] coding for MIMO systems has been an active 
area of research and such codes have been christened Space- 
Time Codes (STC). The primary difference between coded 
modulation (used for SISO, SIMO) and space-time codes is 
that in coded modulation the coding is in time only while in 
space-time codes the coding is in both space and time and 
hence the name. Space-time Codes (STC) can be thought of 
as a signal design problem at the transmitter to realize the 
capacity benefits of MIMO systems [1], [2], though, several 
developments towards STC were presented in [3], [4], [5], [6], 
[7] which combine transmit and receive diversity, much prior 
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to the results on capacity. Formally, a thorough treatment of 
STCs was first presented in [8] in the form of trellis codes 
(Space-Time Trellis Codes (STTC)) along with appropriate 
design and performance criteria. 

The decoding complexity of STTC is exponential in band- 
width efficiency and required diversity order. Starting from 
Alamouti [12], several authors have studied Space-Time Block 
Codes (STBCs) obtained from Orthogonal Designs (ODs) 
and their variations that offer fast decoding (single-symbol 
decoding or double-symbol decoding) over quasi-static fading 
channels [9] -[20], [21]-[27]. But the STBCs from ODs are a 
class of codes that are amenable to single-symbol decoding. 
Due to the importance of single-symbol decodable codes, 
need was felt for rigorous characterization of single-symbol 
decodable linear STBCs. 

Following the spirit of [11], by a linear STBC' we mean 
those covered by the following definition. 

Definition 1 ( Linear STBC): A linear design, S", is a L x 
N matrix whose entries are complex linear combinations of K 
complex indeterminates x^. = a;^/ + jx^g, fc = 0, • • • , — 1 
and their complex conjugates. The STBC obtained by letting 
each indeterminate to take all possible values from a complex 
constellation A is called a linear STBC over A. Notice that 
S is basically a "design"and by the STBC {S, A) we mean 
the STBC obtained using the design S with the indeterminates 
taking values from the signal constellation A. The rate of the 
code/design^ is given by K/L symbols/channel use. Every 
linear design S can be expressed as 

K-l 

5* = ^ XklA2k + XkQA2k+l (1) 
fc=0 

where {Ak\^J^Q^ is a set of complex matrices called weight 
matrices of S. When the signal set A is understood from the 
context or with the understanding that an appropriate signal 
set A will be specified subsequently, we will use the terms 
Design and STBC interchangeably. 

Throughout the paper, we consider only those linear STBCs 
that are obtained from designs. Linear STBCs can be decoded 
using simple linear processing at the receiver with algorithms 
like sphere-decoding [38], [39] which have polynomial com- 
plexity in, A^, the number of transmit antennas. But STBCs 
from ODs stand out because of their amenability to very 
simple (linear complexity in N) decoding. This is because the 
ML metric can be written as a sum of several square terms, 

'Also referred to as a Linear Dispersion code [36] 

-Note that if the signal set is of size 2* the throughput rate R in bits per 
second per Hertz is related to the rate of the design 7?. as i? = TZb. 
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each depending on at-most one variable for OD. However, the 
rates of ODs is restrictive; resulting in search of other codes 
that allow simple decoding similar to ODs. We call such codes 
"single-symbol decodable". Formally 

Definition 2 (Single-symbol Decodable (SD) STBC): A 
Single-symbol Decodable (SD) STBC of rate K/L in K 
complex indeterminates Xk — Xki + jxkQ, k — Q, - ■ ■ ,K — 1 
is a linear STBC such that the ML decoding metric can be 
written as a square of several terms each depending on at 
most one indeterminate. 

Examples of SD STBCs are STBCs from Orthogonal Designs 
of [9]. 

In this paper, we first characterize all linear STBCs that 
admit single-symbol ML decoding, (not necessarily full-rank) 
over quasi-static fading channels, the class of Single-symbol 
Decodable Designs (SDD). Further, we characterize a class of 
full-rank SDDs called Full-Rank SDD (FSDD). 

Fig. n] shows the various classes of SD STBCs identified in 
this paper. Observe that the class of FSDD consists of only 

• an extension of Generalized Linear Complex Orthogonal 
Design (GLCOD^) which we have called Unrestricted 
Full-rank Single-symbol Decodable Designs (UFSDD) 
and 

• a class of non-UFSDDs called Restricted Full-rank 
Single-symbol Decodable Designs (RFSDD)"*. 

The rest of the material of this paper is organized as follows: 
In section|n]the channel model and the design criteria for both 
quasi-static and rapid-fading channels are reviewed. A brief 
presentation of basic, well known results concerning GLCODs 
is given in Section |lll| In Section IIVI we characterize the 
class SDD of all SD (not necessarily full-rank) designs and 
within the class of SDD the class FSDD consisting of full- 
diversity SDD is characterized. Section |V] deals exclusively 
with the maximal rate of square designs and construction of 
such maximal rate designs. 

In Section IVII we generalize the construction of square 
RFSDDs given in Subsection IIV-BI and give a formal defini- 
tion for Co-ordinate Interleaved Orthogonal Designs (CIOD) 
and its generalization. Generalized Co-ordinate Interleaved 
Orthogonal Designs (GCIOD). This generalization is basically 
a construction of RFSDD; both square and non-square and 
results in construction of various high rate RFSDDs. The 
signal set expansion due to co-ordinate interleaving is then 
highlighted and the coding gain of GCIOD is shown to be 
equal to what is defined as the generalized co-ordinate product 
distance (GCPD) for a signal set. A special case of GCPD, 
the co-ordinate product distance (CPD) is derived for lattice 
constellations. We then show that, for lattice constellations, 

'GLCOD is the same as the GeneraUzed Linear Processing Complex 
Orthogonal Design of [9]-the word "Processing" has nothing to be with the 
linear processing operations in the receiver and means basically that the entiies 
are linear combinations of the variables of the design. Since we feel that it 
is better to drop this word to avoid possible confusion we call it GLCOD. 
GLCOD is formally defined in Definition l3l 

^The word "Restricted" reflects the fact that the STBCs obtained from 
these designs can achieve full diversity for those complex constellations that 
satisfy a (trivial) restriction. Likewise, "Unrestricted" reflects the fact that the 
STBCs obtained from these designs achieve full diversity for all complex 
constellations. 



GCIODs have higher coding gain as compared to GLCODs. 
Simulation results are also included for completeness. The 
maximum mutual information (MMI) of GCIODs is then 
derived and compared with that of GLCODs to show that, 
except for iV = 2, CIODs have higher MMI. In short, this 
section shows that, except for N = 2 (the Alamouti code), 
CIODs are better than GLCODs in terms of rate, coding gain 
and MMI. 

In section IVIII we study STBCs for use in rapid-fading 
channels by giving a matrix representation of the multi-antenna 
rapid-fading channels. The emphasis is on finding STBCs that 
allow single-symbol decoding for both quasi-static and rapid- 
fading channels as BER performance such STBCs will be 
invariant to any channel variations. Therefore, we characterize 
all linear STBCs that allow single-symbol ML decoding when 
used in rapid-fading channels. Then, among these we identify 
those with full-diversity, i.e., those with diversity L when the 
STBC is of size L x N,{L > N), where is the number 
of transmit antennas and L is the length of the code. The 
maximum rate for such a full-diversity, SD code is shown to 
be 2 /L from which it follows that rate-one is possible only for 
2 Tx. antennas. The co-ordinate interleaved orthogonal design 
(CIOD) for 2 Tx (introduced in Section llVt is shown to be 
one such rate-one, full-diversity and SD code. (It turns out that 
Alamouti code is not SD for rapid-fading channels.) Finally, 
Section lVlIII consists of some concluding remarks and a couple 
of directions for further research. 

II. Channel Model 

In this section we present the channel model and review the 
design criteria for both quasi-static and rapid-fading channels. 
Let the number of transmit antennas be N and the number of 
receive antennas be M. At each time slot t, complex signal 
points. Sit, « = 0, 1, • • • , — 1 are transmitted from the N 
transmit antennas simultaneously. Let hijt = aijte^^'^* denote 
the path gain from the transmit antenna i to the receive antenna 
j at time t, where j = ^/—l. The received signal Vjt at the 
antenna j at time t, is given by 



N-l 
i=0 



(2) 



j = 0, • • • , M — 1; t — 0, - ■ ■ ,L — 1. Assuming that perfect 
channel state information (CSI) is available at the receiver, the 
decision rule for ML decoding is to minimize the metric 



L-l M-1 

EE 

t=0 1=0 



ijt ^it 



(3) 



over all codewords. This results in exponential decoding 
complexity, because of the joint decision on all the symbols sa 
in the matrix S. If the throughput rate of such a scheme is R in 
bits/sec/Hz, then 2^^ metric calculations are required; one for 
each possible transmission matrix S. Even for modest antenna 
configurations and rates this could be very large resulting in 
search for codes that admit a simple decoding while providing 
full diversity gain. 
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A. Quasi-Static Fading Channels 

For quasi-static fading channels hijt = hij and (|3 can be 
written in matrix notation as, 

JV-l 

Vjt = ^ hijSit+rijt, j = 0, •••,M-1; < = 0, •••,L-1. 

i=0 

(4) 

In matrix notation, 

V = SH + W (5) 

where V e £Lxm denotes the complex field) is the 
received signal matrix, S G C^^^ is the transmission matrix 
(codeword matrix), H G £^nxm denotes the channel matrix 
and W S £^LxM jj^g entries that are Gaussian distributed 
with zero mean and unit variance and also are temporally and 
spatially white. In V, S and W time runs vertically and space 
runs horizontally. The channel matrix H and the transmitted 
codeword S are assumed to have unit variance entries. The 
ML metric can then be written as 

M(S) = tr ((V - SH)^(V - SH)) . (6) 

This ML metric (|SJl results in exponential decoding complexity 
with the rate of transmission in bits/sec/Hz. 

1 ) Design Criteria for STC over quasi-static fading chan- 
nels: The design criteria for STC over quasi-static fading 
channels are [8]: 

• Rank Criterion: In order to achieve diversity of rM, the 
matrix B(S, S) = S — S has to be full rank for any two 
distinct codewords S, S. If B(S, S) has rank N, then the 
STC achieves full-diversity. 

• Determinant Criterion: After ensuring full diversity the 
next criteria is to maximize the coding gain given by, 

A(S,S) = min|(S-S)(S-S)^iy'' (7) 
s,s 

where \A\^ represents the product of the non-zero eigen 
values of the matrix A. 

2) Design Criteria for STC over Rapid-Fading Channels:: 
We recall that the design criteria for rapid-fading channels are 
[8]: 

• The Distance Criterion : In order to achieve the diversity 
rM in rapid-fading channels, for any two distinct code- 
word matrices S and S, the strings sot, sit, ■ ■ • , S(Ar-i)t 
and Sot, su, • • • , S(jv-i)t must differ at least for r values 
ofO < t < L — 1. (Essentially, the distance criterion 
implies that if a codeword is viewed as a L length vector 
with each row of the transmission matrix viewed as a 
single element of C^, then the diversity gain is equal to 
the Hamming distance of this L length codeword over 
C^. 

• The Product Criterion : Let V(S, S) be the indices of the 
non-zero rows of S — S and letjst — Stp — ^fj^^ \sit — 
Sitp, where St is the t-th row of S, < i < L — 1. Then 
the coding gain is given by 

min TT |st-stp. 

tev(s.s) 

The product criterion is to maximize the coding gain. 



III. Generalized Linear Complex Orthogonal 
Designs (GLCOD) 

The class of GLCOD was first discovered and studied in 
the context of single-symbol decodable designs by coding 
theorists in [9], [11], [19], [17], [51]. It is therefore proper 
to recollect the main results concerning GLCODs before 
the characterization of SSD. In this section we review the 
definition of GLCOD and summarize important results on 
square as well as non-square GLCODs from [9], [11], [19], 
[17], [51]. 

Definition 3 (GLCOD): A Generalized Linear Complex 
Orthogonal Design (GLCOD) in k complex indeterminates 
xi,X2,--- ,Xk of size N and rate TZ = k/p, p > N is a 
p X N matrix Q, such that 

• the entries of are complex linear combinations of 
0, ±Xi, i — 1, - ■ ■ ,k and their conjugates. 

• 8^0 — D, where D is a diagonal matrix whose entries 
are a linear combination of i = 1, • ■ ■ i ^ with all 
strictly positive real coefficients. 

If k=N=p then Q is called a Linear Complex Orthogonal 
Design (LCOD). Furthermore, when the entries are only from 
{0, ±xi, ±X2, • • • , ±a;/c}, their conjugates and multiples of j 
then Q is called a Complex Orthogonal Design (COD). STBCs 
from ODs are obtained by replacing Xi by Si and allowing Si 
to take all values from a signal set A. A GLCOD is said to 
be of minimal-delay if N — p. 

Actually, according to [9] it is required that D = 
J2'i=i which is a special case of the requirement 

that D is a diagonal matrix with the conditions in the 
above definition. In other words, we have presented a 
generalized version of the definition of GLCOD of [9]. Also 
we say that a GLCOD satisfies Equal- Weights condition if 

The Alamouti scheme [12], which is of minimal-delay, full- 
rank and rate-one is basically the STBC arising from the size 
2 COD. 

Consider a square GLCOD^, S = X^s^o^ ^ki ^2fe + 
XkQA2k+i- The weight matrices satisfy, 

Al^Ak^Vk, fc = 0,---,2X-l (8) 
Aj^Ak + A^Ai =0, 0<k^l<2K ^1. (9) 

where 2?^ is a diagonal matrix of full-rank for all k. Define 

^ 1/2 

Bk = AkV^ . Then the matrices Bk satisfy (using the 
results shown in [51]) 

Bj^Bk^lN, k ^ {),■■■ ,2K-\ (10) 
B^Bk + B^Bi^Q, Q<k^l<2K-l (11) 

and again defining 

Ck = B^Bk, k^Q,--- ,2K-\, (12) 

we end up with Co = In and 

C^ = -Cfc, fc = l,---,2X-l (13) 
Cl^Ck + C^Ci =0, l<k^l<2K-l. (14) 

rate-1, square GLCOD is referred to as complex linear processing 
orthogonal design (CLPOD) in [9]. 
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The above normalized set of matrices {Ci, • • • , C2K-1} con- 
stitute a Hunvitz family of order N [28]. Let H {N) — l denote 
the number of matrices in a Hurwitz family of order N, then 
the Hurwitz Theorem can be stated as 

Theorem 1 (Hurwitz [28]): If 2''b, b odd and a,b> 
then 

H{N) <2a + 2. 
Observe that H (TV) = 2K. An immediate consequence of the 
Hurwitz Theorem are the following results: 

Theorem 2 (Tarokh, Jafarkhani and Calderbank [9]): 
A square GLCOD of rate-1 exists iff N = 2. 

Theorem 3 (Trikkonen and Hottinen [11]): The 
maximal rate, 7?, of a square GLCOD of size TV = 2''6, 6 odd, 
satisfying equal weight condition is 



n 



1 



N 



This result was generalized to all square GLCODs in [51] 
using the theorem: 

Theorem 4 (Khan and Rajan [51]): With the Equal- 
Weights condition removed from the definition of GLCODs, 
an N X N square (GLCOD), £c in variables xq,- 
exists iff there exists a GLCOD Cr such that 



,XK-1 



)I. 



(15) 



C^Ca = i\xo\' + --- + \XK-1\ 

Hence we have the following corollary. 

Corollary 5 (Khan and Rajan [51]): Let N = 2"6 where 
b is an odd integer and a = 4c + d, where < d < c and 
c > 0. The maximal rate of size N, square GLROD without 
the Equal- Weights condition satisfied is ^^jr- and of size N, 



square GLCOD without the Equal-Weights condition satisfied 
is 

An intuitive and simple realization of such GLCODs based on 
Josefiak's realization of the Hurwitz family, was presented in 
[19] as 

Construction 3.1 (Su and Xia [19]): Let Gi{xq) = XqIi, 
then the GLCOD of size 2^, G2k{xo,xi, ■ ■ ■ ,xk), can be 
constructed iteratively for K = 1,2,3, ■■ ■ as 



G2k{xo,Xi, ■ ■ ■ ,Xk) = 

G2K-l{xQ,Xl,- ■ ■ ,XK-l) XKh^-i 

-Xxhii-i Gl^K-i{xo,xi,- ■ ■ ,XK-i) 

While square GLCODs have been completely characterized 
non-square GLCODs are not well understood. The main results 
for non-square GLCODs are due to Liang and Xia. The 
primary result is 

Theorem 6 (Liang and Xia [16]): A rate 1 GLCOD ex- 
ists iff iV 2. 

This was further, improved later to. 

Theorem 7 (Su and Xia [19]): The maximum rate of 
GCOD (without linear processing) is upper bounded by 3/4. 
Xue bin-Liang [17] gave the construction of maximal rates 
GCOD 

Theorem 8 (Liang [17]): The maximal rate of a GCOD 
for g N transmit antennas is given hy TZ = ^^^^ where 
m = lN/2\ . 

The maximal rate and the construction of such maximal rate 
non-square GLCODs for N > 2 remains an open problem. 



IV. Single-symbol Decodable Designs 

In the first part of this section we characterize all STBCs 
that allow single-symbol ML decoding in quasi-static fading 
channel and using this characterization define Single-symbol 
Decodable Designs (SDD) in terms of the weight matrices and 
discuss several examples of such designs. In the second part, 
we characterize the class FSDD and classify the same. 



A. Characterization of SD STBCs 

Consider the matrix channel model for quasi-static fading 
channel given in ^ and the corresponding ML decoding 
metric (|6|l. For a linear STBC with K variables, we are 
concerned about those STBCs for which the ML metric (|6j can 
be written as sum of several terms with each term involving 
at-most one variable only and hence SD. 

The following theorem characterizes all linear STBCs, in 
terms of the weight matrices, that will allow single-symbol 
decoding. 

Theorem 9: For a Unear STBC in K variables, S — 



J2k=o ^kiA2k + XkQA2k+i, the ML metric, M{S) defined 
in (|6j decomposes as AI{S) — ^^^q Mk{xk) + Mc where 
Ale — —{K — l)tr (y^y) is independent of all the variables 
and Mi.{xi.) is a function only of the variable Xk, iff* 



\fl^k,k- 



1 if k is even 
1 if k is odd 



(17) 



Proof From (|6} we have 
M{S) = tr (V^V) - tr ((5H)^V) 



- tr 



(V^5H)) 



Observe that tr (V^V) is independent of S. The next two 
terms in M{S) are functions of S, S'^ and hence hnear in 
Xki, XkQ. In the last term. 



S^S = 



(16) 



K-l 

k=0 

A'-l A'-l 

+ E E 

k=0 l=k+l 



{A^kA2kxti + A^k^^A2k+ixiQ) 



{A^,A2i + A^iA2k 



K-1 K-1 

(^2fc^ 



E 



fc=0 l=k+l 
K-1 K-1 



l^2fc^ 



f^kQXlQ 



E E(^?fc^2;+i 



A^i+iA2k)xkixiQ. (18) 



fc=0 (=0 



(a) Proof for the "if part": If is satisfied then (I18> reduces 

'The condition 1171 can also be given as 



due to the 
tr{(V- SH)(V- 



-A,Al 







^l^k,k- 



1 if k is even 
1 if k is odd 



identity tr { (V - SH)'" (V - SH 
SH)^ } when 5" is a square matrix 
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to 



K-1 



S^S = i^2kA2kxl + A^,^,A2k+ixlQ 

+ (^?fc^2fe+l + ^?fe+1^2fe) XklXkq) 

= J2 ^kTk, where 

fc=0 

Tk = A2kXkI + A2k+lXkQ 



(19) 
(20) 



and using linearity of the trace operator, M{S) can be written 
as 



M{S) 



tr(V«V)-^{tr((T,H)«V) 



fe=0 



-tr (V^TfeH) + tr {Tk^nUn^)} 

= ^ ||V - {A^kXkl + A2k+iXkQ)-af +H^l) 

where = -{K - iy'^V) and ||.|| denotes the 

Frobenius norm. 

(b) Proof for the "only if part": If (I17> is not satisfied for 

any Au^ , Ai^ ^ki^h then 

M{S) = ^||V-(A2fca:fc/+A2fc+ia;fcQ)H||' 

k 

+tr A,, + A^AfcjH'^H) y + AU22) 

where 

{2;(fci/2)/a;(;i/2)/ if both ki,li are even 

a;((fci-i)/2)ga;((ii_i)/2)Q if both fci,Zi are odd 
a;((fci-i)/2)Qa;(ii/2)/ if fci odd, ?i even. 

Now, from the above it is clear that M{S) can not be 
decomposed into terms involving only one variable. ■ 

It is important to observe that Ml\ implies that it is not 
necessary for the weight matrices associated with the in-phase 
and quadrature-phase of a single variable (say fc-th) to satisfy 
the condition Al^j^^^A2k + Al^^A2k+i = 0. Since Al^f,^^A2k + 
^2fe^2fc+i is indeed the coefficient of XkiXkQ in S'^S', this 
implies that terms of the form XkiXkQ can appear in S'^S' 
without violating single-symbol decodability. An example of 
such a STBC is given in Example 14. II 

Example 4.1: Consider 



X0Q+}XiQ Xoz+jxi/ 



The corresponding weight matrices are given by 

A2 



(23) 



" 1 





,^1 = 


" 


1 





1 


1 





j 




" 

j 




" 

_ j 


j 





A^j^A2k+i 7^ for fc = as well as fc = 1. Explicitly, 

A^A^+A^Ao ^ (24) 
A?A3 + A^A2 ^ 



A^A2 + A^Ao = 
^^A^ + A^Ao = 

A^A, - 



A^A2 

A^A2 



^3 



(25) 
(26) 
(27) 

(28) 

12+/12/11 - 0. (29) 
Remark 10: However note that for the SD STBC in Exam- 
ple O 

det {(5 - 5)^(5 -S)] = [{Ax^i - A.toq)' + 

(Axi/ - Axiq)2] [(Axo7 + Axoq)2 + (Axi7 + Axiq)2] 

where Xi — Xi — ISxu +]lSxiQ. If we set A.ti/ = Axig = 
we have 



dot [{S- S)'^{S - S)] = [(A2a;o/ - A^xoq)' 



(30) 



which is maximized (without rotation of the signal set) when 
either A^xq/ = or A^xqq = 0, i.e. the fc-th indeterminate 
should take values from a constellation that is parallel to 
the "real axis" or the "imaginary axis". Such codes are 
closely related to Quasi-Orthogonal Designs (QOD) and the 
maximization of the corresponding coding gain with signal set 
rotation has been considered in [58], [59]. 

Henceforth, we consider only those STBCs S = 
Yjk=o XkiA2k + XkQA2k+i, which have the property that 
the weight matrices of the in-phase and quadrature com- 
ponents of any variable are orthogonal, that is 



A^ 



0<fc</v-l (31) 



i2fc^2A;+l -I- ^2fe+1^2fe - 0, 

since all known STBCs satisfy ^31) and we are able to 
tract and obtain several results concerning fuU-rankness, 
coding gain and existence results with this restriction. 

Theorem |9] for this case specializes to: 

Theorem 11: For a linear STBC in K complex variables, 
S = Yl,k=a XkiA2k + XkQA2k+i satisfying the necessary 
condition ^^^,^2^+1 + ^^^+1^2^ = 0, < fc < i^T - 1, the 
ML metric, M{S) defined in (|6j decomposes as M{S) = 
Ek=o Mk{xk) + M, where M, = -{K - l)tr [V^V), iff 

A^A,+A^Afe=0, < fc 7^ ; < 2is:- 1. (32) 
We also have 

Proposition 12: For a linear STBC in K complex variables, 
S = Yl,k=o ^ki A2k + XkQA2k+i satisfying the necessary 
condition A^,,A2k+i+ A^^+i^'^k = 0, < fc < K-l, the ML 
metric, M{S) defined in (|6} decomposes as M{S) = ^ 
Mk{xk) + Mc where = -{K - l)tr (V'^V), iff 

tr (AfcHH^A^ + AiHli^Al^) =0, G<fc7^/<2if-l. 

(33) 

If, in addition, S is square (N = L), then (13 3> is satisfied if 
and only if 

AuA\^ + AiA^ 0<k^l<2K -1. (34) 

Proof: Using the identity. 



k=0 



and it is easily verified that O is satisfied and ^^^+1^2^ + tr (^(V - SH)'"(V - SH)) = tr ((V - SH)(V - SH 
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can be written as 

M(5) = tr (VV") - tr ((S'H)V^) - tr (V(S'H)^) 
+tr (5HH^5") . 

Observe that tr(VV^) is independent of S. The next two 
terms in M{S) are functions of S, 5^ and hence linear in 
Xki, XkQ- In the last term, 



kQ) 



K-l 
fe=0 

k=0 l=k+l 
K-1 K-1 

k=0 l=k+l 
K-1 K-1 

+ (^2feB^+i + B2i+iBli)xkiXiQ (35) 

k=0 1=0 

where Bk = AkH (a) Proof for the "if part": If (jSSjl is satisfied 
then i35\ reduces to 

K-l 



which is satisfied iff AkAj^ + AiA]^ = 0, Q < k ^ I < 

Examples of SD STBCs are those from OD, in-particular 
the Alamouti code. The following example gives two STBCs 
that are not obtainable as STBCs from ODs. 

Example 4.2: For N ~ K — 2 consider 

xoi+jxiQ 

XlI+jXQQ 

The corresponding weight matrices are given by 





(40) 



Ao 



A. 








1 



,Ai 



,A, 



J 






Similarly, for N = K = 4 consider the design given in J41t . 
The corresponding weight matrices are 



S^S = J2 (^2fcHH«A?fea;L + A^k+iIlH^ A^,^, 

k=0 
K-l 

= TkT'^, where 



k=0 



(36) 

Tk = {A2kXki + A2k+iXkQ)ll (37) 

and using linearity of the trace operator, M (S) can be written 

as 

K-l 

M{S) = tr (VV'") - {tr (^feV^) - tr (VTfc") 



Ao 



A2 - 



fc=0 



A. 



+tr (TfcT, 

= V ||V - {A2kXki + A2k+iXkQ)nf +Af^38) 

where Mc = -{K - l)ti: (V'^V) and ||.|| denotes the 
Frobenius norm. 

(b) Proof for the "only if part": If (I33t is not satisfied for 

any Ak^ , Ai^ .ki^h then 

M{S) = ^||V-(A2feXfe/ + A2fe+iXfeQ)H||' 



1 










" 





1 



































_ 







1 








-1 




















































' 
























1 
















1 









































1 










-1 






,^1 



,^3 = 










j 




j 

-j 
















.A, 



where 



X(ki/2)iX(ij2)i if both ki,li are even 

X((ki-i)/2)QX((i^-i)/2)Q if both fci, ^1 are odd 
X((ki-i)/2)QX(i^/2)i if fci odd, h even. 

Now, from the above it is clear that M{S) can not be 
decomposed into terms involving only one variable. 
For square 5, ( I33> can be written as 

tr(HH" {AfeA]^ + AjA^}) = 0, Q<kT^l<2K-\ 

(39) 



It is easily seen that the two codes of the above example 
are not covered by GLCODs and satisfy the requirements of 
Theorem ll H and hence are SD. These two STBCs are instances 
of the so called Co-ordinate Interleaved Orthogonal Designs 
(CIOD), which is discussed in detail in Section IVII and a 
formal definition of which is Definition These codes apart 
from being SD can give STBCs with full-rank also when the 
indeterminates take values from appropriate signal sets- an 
aspect which is discussed in detail in Subsection II V-BI and in 
Section ED 

B. Full-rank SDD 

In this subsection we identify all full-rank designs with in 
the class of SDD that satisfy ( I32> . calling them the class of 
Full-rank Single-symbol Decodable Designs (FSDD), charac- 
terize the class of FSDD and classify the same. Towards this 
end, we have for square (N = L) SDD 

Proposition 13: A square SDD S — XkiA2k + 

XkQA2k+i, exists if and only if there exists a square SDD, 



S ■ 



A 



XkiA 



2k 



, Ai + A^Ak 



+ XkQA2k 

0, k^l, 



1 such that 
and 



A^Ak 
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-Xil + jX3Q 






Xll +JX3Q 









X2I + ja:oQ 



X2I 



(41) 



where Vk is a diagonal matrix. 

Proof: Using (|32} and ( I34l i repeatedly we get 

A]^AuAYAi - A^{~AiA'^)Ai^{A^Au)A^Ai 
= AYAki-A^Ak) - 

which implies that the set of matrices {-A^ AkY^J^^^ forms 
a commuting family of Hermitian matrices and hence can be 
simultaneously diagonalized by a unitary matrix, U . Define 
Ak = AkU"^, then S = Y.k=o^ XkiA2k+XkQA2k+i is a linear 
STBC such that = 0,Vfc ^ l,A]^Ak = Vu^h 

where Vk is a diagonal matrix. For the converse, given S, 
S — SU where [/ is a unitary matrix. ■ 
Therefore for square SDD, we may, without any loss of 
generality, assume that S'^S is diagonal. To characterize non- 
square SDD, we use the following 

Property 4.1 (Observation 7.1.3 of [65]): Any non- 
negative linear combination of positive semi-definite matrices 
is positive semi-definite. 
Propertv l4. 1 1 when applied to a SDD yields 

Property 4.2: For a SDD, S = Y.k=o XkiA2k + 
XkQA2k+i, the matrix 5^5* is positive semi-definite and 
A^Ak, Vfc are positive semi-definite. 

Using proDertv l4.2l we have the following necessary condition 
for a SDD to have full-diversity. 

Proposition 14: If an SDD, S = Y^k=o XkiA2k + 
XkQA2k+i, whose weight matrices Ak satisfy 

Aj^Ai+A^Ak^O, yk^l (42) 

achieves full-diversity then AJ^^A2k+ ^k+i^'2-k+i is full-rank 
for all k = 0, 1, • • • ,K — 1. In addition if S is square then 
the requirement specializes to I?2fc + 2^2fe+i being full-rank 
for all fc = 0, 1, • • • , X — 1, where the diagonal matrices T>i 
are those given in Proposition 1 131 

Proof: The proof is by contradiction and in two parts 
corresponding to whether 5 is square or non-square. 
Part 1: Let S* be a square SDD then by Proposition^] without 
loss of generality, Ak — Vk^k. Suppose V2k + T^2k+i, 
for some he [0, iiT - 1], is not full-rank. Then S'^S = 
J2k=o '^2kxlj + T>2k+ixlQ. Now for any two transmission 
matrices S, S that differ only in Xk, the difference matrix 
B{S, S) = S-S, will not be full-rank as B'^{S, S)B{S, S) = 
T>2k{xki - Xki)^ + T>2k+i{xkQ - XkqY is not full-rank. 
Part 2: The proof for non-square SDD, S, is similar to 
the above except that B'^{S,S) B{S,S) = A'^^,A2k{xki - 
XkiY + Ay-f^^^A2k+i{xkQ - XkqY where A]^ Ak are positive 
semi-definite. Since a non-negative linear combination of 
positive semi-definite matrices is positive semi-definite, for 
full-diversity it is necessary that A^A2fe + Ayj^_^_^A2k+i is 
full-rank for all fc = 0, 1, • • • , - 1. ■ 
Towards obtaining a sufficient condition for full-diversity, 
we first introduce 



Definition 4 (Co-ordinate Product Distance (CPD)): The 
Co-ordinate Product Distance (CPD) between any two signal 
points u = ui + juQ and v = vj + jvq, u ^ v, in the signal 
set A is defined as 

CPD{u,v)^\ui-vi\\uq-vq\ (43) 

and the minimum of this value among all possible pairs is 
defined as the CPD of A. 

Remark 15: The idea of rotating QAM constellation was 
first presented in [60] and the term "co-ordinate interleaving" 
as also "Co-ordinate Product Distance" was first introduced by 
Jelicic and Roy in [42], [43] in the context of TCM for fading 
channels. This concept of rotation of QAM constellation was 
extended to multi-dimensional QAM constellations in [61], 
[62] at the cost of the decoding complexity. However, for 
the two-dimensional case there is no increase in the decoding 
complexity as shown in [40], [41]. 

Theorem 16: A SSD, S = Y.k~^ XkiA2k+XkQA2k+i 
where Xk take values from a signal set A,yk, satisfying the 
necessary condition of Proposition achieves full-diversity 
iff 

(i) either A^Ak is of full-rank for all k or (ii) the 
CPD of ^ 0. 

Proof: Let be a square SDD satisfying the 
necessary condition given in Theorem We have 

B^{Sj)BiSj)=j:^^^V2k+i ixki'Xkif+V2k+iixkQ- 
XkQY ■ Observe that under both these conditions the difference 
matrix B{S, S) is full-rank for any two distinct S, S. Con- 
versely, if the above conditions are not satisfied then for exist 
distinct S, S such that -8(5', S) is not full-rank. The proof is 
similar when S" is a non-square design. ■ 
Examples of FSDD are the GLCODs and the STBCs of 
Example 14.21 

Note that the sufficient condition (i) of Theorem ^] is 
an additional condition on the weight matrices whereas the 
sufficient condition (ii) is a restriction on the signal set A and 
not on the weight matrices Ak- Also, notice that the FSDD that 
satisfy the sufficient condition (i) are precisely an extension of 
GLCODs; GLCODs have an additional constraint that Aj^Ak 
be diagonal. 

An important consequence of Theorem^] is that there can 
exist designs that are not covered by GLCODs offering full- 
diversity and single-symbol decoding provided the associated 
signal set has non-zero CPD. It is important to note tliat 
wlienever we have a signal set with CPD equal to zero, by 
appropriately rotating it we can end with a signal set with 
non-zero CPD. Indeed, only for a finite set of angles of 
rotation we will again end up with CPD equal to zero. So, 
the requirement of non-zero CPD for a signal set is not at 
all restrictive in real sense. In Section fvTl we find optimum 
angle(s) of rotation for lattice constellations that maximize the 
CPD. 
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For the case of square designs of size N with rate-one it is 
shown in Section IVl that FSDD exist for N — 2,4 and these 
are precisely the STBCs of Example 14.21 and the Alamouti 
code. 

For a SDD, when Aj^Ak is full-rank for all k, corresponding 
to Theorem^]with the condition (i) for full-diversity satisfied, 
we have an extension of GLCOD in the sense that the STBC 
obtained by using the design with any complex signal set 
for the indeterminates results in a FSDD. That is, there is no 
restriction on the complex signal set that can be used with 
such designs. So, we define. 

Definition 5 (Unrestricted FSDD (UFSDD)): A FSDD is 
called an Unrestricted Full-rank Single-symbol Decodable De- 
sign (UFSDD) if A^Ak is of full-rank for all fc = 0, • • • , 2K~ 
1. 

Remark 17: Observe that for a square UFSDD S, Aj^Ak = 
Vk is diagonal and hence UFSDD reduces to square GLCOD. 
For non-square designs, GLCOD is a subset of UFSDD. 
Also the above extension of the definition of GLCODs was 
hinted in [19] where they observe that A^Ak can be positive 
definite. However it is clear from our characterization that 
such a generalization does not result in any gain for square 
designs. For non-square designs existence of UFSDDs that are 
not GLCODs or unitarily equivalent to GLCODs is an open 
problem. 

The FSDD that are not UFSDDs are such that Al^^,A2k and/or 
^2l-+i^2fe+i is not full-rank for at least one k. (The CIOD 
codes of Example 14.21 are such that I?2fc +2^2fc+i is full-rank 
Vfc and Vk is not full-rank for all k.) We call such FSDD 
codes Restricted Full-rank Single-symbol Decodable Designs 
(RFSDD), since any full-rank design within this class can be 
there only with a restriction on the complex constellation from 
which the indeterminates take values, the restriction being that 
the CPD of the signal set should not be zero. Formally, 

Definition 6 (Restricted FSDD (RFSDD)): A Restricted 
Full-rank Single-symbol Decodable Designs (RFSDD) is a 
FSDD such that A^Ak is not full-rank for at least one k 
where fc = 0, • • • , 2K — 1 and the signal set, from which the 
indeterminates take values from, has non-zero CPD. 
Observe that the CIODs are a subset of RFSDD. Figure [T] 
shows all the classes discussed so far, viz., SDD, FSDD, 
RFSDD, UFSDD. In Section |Vl we focus on the square 
RFSDDs as square UFSDD have been discussed in Section 

Iml 

V. Existence of Square RFSDDs 

The main result in this section is that there exists square 
RFSDDs with the maximal rate || for = 2° antennas 
whereas only rates up to ^tl is possible with square 
GLCODs with the same number of antennas. The other 
results are: (i) rate-one square RFSDD of size N exist, iff 

= 2, 4 and (ii) a construction of RFSDDs with maximum 
rate from GLCODs. 

Let S = Y.k=o ^kiA2k + XkQA2k+i be a square RFSDD. 
We have, 

Aj^Ak=Vk, k^0,---,2K-l (44) 
A^Ak + Aj^Ai = 0, 0<k^l<2K-l (45) 



where Vk, k — 0, - ■ ■ , 2K ~1 are diagonal matrices with non- 
negative entries such that I?2fe + T^2k+i is full-rank Vfc. First 
we show that for a rate-one RFSDD, 2, 4 or 8. 

Theorem 18: If S* is a size N square RFSDD of rate-one, 
then iV = 2, 4 or 8. 

Proof: Let Bk = A2k + A2k+i, A: = 0, • • • , A' - 1, 

then 

B'^Bk='Dk=V2k+V2k+i, k^O,---,K-l (46) 
Bl^Bk + BJ^Bi =0, 0<k^l<K-l. (47) 

Observe that T>k is of full-rank for all k. Define Ck — 
BkVf^ ' . Then the matrices C'k satisfy 

Cj^Ck = lN: k = Q,---,K-\ (48) 
Cl^Ck + C^Ci =0, Q<k^l<K~l. (49) 

Define 

Ck^C^Ck, k = 0,---,K-l, (50) 

then Co — In and 

C'^ = ~Ck, k^l,---,K~l (51) 
CPCk + Cj^Ci =0, l<k^l<K -1. (52) 

The normalized set of matrices {Ci, • • • , Ck-i} constitute a 
Hurwitz family of order N [28] and for N = 2"6, b odd and 
a,b > the number of such matrices if — 1 is bounded by 
[28] 

K <2a + 2. 

For rate-one, RFSDD (K = N), the inequality can be satisfied 
only for TV = 2,4 or 8. ■ 
Therefore the search for rate-one, square RFSDDs can be 
restricted to = 2,4,8. The rate 1, RFSDDs for A^ = 2,4 
have been presented in Example 14.21 We will now prove that 
a rate-one, square RFSDD for A^ = 8 does not exist. Towards 
this end we first derive the maximal rates of square RFSDDs. 

Theorem 19: The maximal rate, TZ, achievable by a square 
RFSDD with A^ = 2"6,6 odd (where a,b > 0) transmit 
antennas is 

2a 

^ = TTT (53) 

K-l 

Proof Let 5" = Y.k=o XkiA2k+XkQA2k+i be a square 
RFSDD. Define the RFSDD 

S'^V] xki C^A2kVo^'^ -^XkQ C^A2k+i'D-^'^ 
— ^ ^ ^ ^ 

fc=o 7; "7' 

where Ck and Vk are defined in the proof of the previous 
theorem. Then the set of matrices {C'j. = A^k + ^2fe+i} i^ 
such that Cq = In and {C^, A: = 1, • • • , if — 1} is a family 
of matrices of order A^ such that 

C'k ^Vo^Vk, l<k<K-l, (54) 
C'l^C'^ + C'f^C'i^Q, 0<k^l<K-l, (55) 

where V^^Vk is diagonal and full-rank for all k. Then we 
have 

A'o + A[ = C'o^In. (56) 
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It is easily that the set of matrices {A'l^} satisfy J44> and 
( I45> . Also, at least one A'/, is not full-rank. Without loss of 
generality we assume that A'q is of rank r < (if this not 
so then exchange the indeterminates and/or the in-phase and 
quadrature components so that this is satisfied). As A'q is of 
rank r, due to ( I44> . n — r columns of Aq are zero vectors. 
Assume that first r columns of A'q are non-zero (If this is not 
the case, we can always multiply all the weight matrices with 
a Permutation matrix such that A'q is of this form) i.e. 



A'q=[B'q ] 



(57) 



where B'q e C^^''. Applying (|45} to A'q and A'^^ and using 
from (I56> and ( I57> . we have 



A'q''{In-A'q) + (In~A'q'')A'q^Q 
2V'q 

[B'q ]^2V'q 



^0 



^q^2V'q 



^0 







(58) 
(59) 

(60) 
(61) 



where B'n is a r x r matrix and full-rank and A'j!^A'i. = 
V'l^, k — 0, - ■ ■ , 2K — 1. Therefore the matrices A'q, A'l are 
of the form 



A' — 



B'l, 




Ir — B'li 




IN- 



(62) 



Let 



Di = 



Dn Di2 

^21 ^22 



be a matrix such that 

A'^Di + D^A'^ = 0, 



0,1 



(63) 



where Dn e C'■^^ £'22 e 
structure of An we have 



-'N—rxN—r 



Substituting the 



A'^D^+D^A'q = 

B'^Dn+D]{B'^^ B'^D^2 
D^2B[i 



0. 



(64) 
(65) 



As B'^^ is full-rank it follows that D12 
structure of A'^ we have 



{Ir - B'^)Dii + D'^^ilr - B'n) 

D21 



-- 0. Substituting the 

^21 



D2 



D^2 



(66) 
D21 - 0. (67) 



It follows that Di is block diagonal and consequently all the 
A'j., 2 < k < 2K — 1 are block diagonal of the form Di as 
they satisfy j63> . Consequently, = A' 2k + A'2k+l^ k = 
1, ■ ■ ■ , K — 1 are also block diagonal of the form = 



Also, from 



where C'^^ £ C 
, (l66t we have 



•<N—rxN—r 



DV„D2 



k2 



^22- 



(68) 



Now, in addition to this block diagonal structure the matrices 
A'^, 2 < k < K — 1 have to satisfy ( I45> among themselves. 



It follows that the two sets of square matrices {C'f.^,k 
0, . . . , A' - 1} and {C'f.^, fc = . . . , AT - 1} satisfy 



^r2 _ 

^ki — " 

^ki^H 



,K-l, 



1,2; 



l<k^l<K-l, 



(69) 



where —Vki are diagonal and full-rank Vfc, i. Define 

(^fc,: = C;,2?,//^ k^l,---,K-l, z = l,2; (71) 
then from Theorem |4] 

k = l,--- ,K-1, i = l,2; (72) 
Ck^Cl, = -Cl,Ck^, l<k^l<K-l, t = ICSS) 

and the sets of square matrices {Cki, k — 1,2, . . . , K — 1} 
and {Ck2, k — 1,2, . . . , K — 1} constitute Hurwitz families 
of order r,N ~ r corresponding to i = 1,2 respectively. Let 
H (N) — 1 be the maximum number of matrices in a Hurwitz 
family of order N, then from the Hurwitz Theorem [28] , 
N = 2^1,1 odd and 

H (N) = 2a + 2. (74) 

Observe that due to the block diagonal structure of C^, K = 
min{_ff (rj , H {N — r^)}. Following the Hurwitz Theorem it 
is sufficient to consider both r,N — r to be of the form 2", 
say 2°i , 2"^ respectively. It follows that K is maximized iff 
r = N-r = 2°-'^N = 2°^'+^. It follows that the maximum 
rate of RFSDD of size = 2" (a = a' + 1) is 



^ 2a 

2" 



(75) 



An important observation regarding square RFSDDs is sum- 
marized in the following Corollary: 

Corollary 20: A maximal rate square RFSDD, S = 
Y.k=Q XkiA2k + XkQA2k+i exists iff both X'2fc and X'2fc+i 
are not full-rank for all k. 

Proof: Immediate from the proof of above theorem. ■ 
An immediate consequence of this characterization of maximal 
rate RFSDDs is: 

Theorem 21: A square RFSDD of rate-one, exists iff N = 
2,4. 

Proof: From (|75} 7^ = 1 iff TV = 2, 4 ■ 

It follows that 

Theorem 22: The maximal rate, TZ, achievable by a square 
FSDD with N ~ 2°-b, b odd (where a,b > 0) transmit antennas 
is 



7^ 



2a 
2°6' 



(76) 



Furthermore square GLCODs are not maximal rate FSDD 
except for N — 2. 

Next we give a construction of square RFSDD that achieves 
the maximal rates obtained in Theorem 1 191 

Theorem 23: A squai-e RFSDD S, of size N, in variables 
Xi,i = Q, ■ ■ ■ ,K—1 achieving the rate of Theorem^] is given 
by 



e(io, • • • ,xk/2) 

'd{xK/2,-- ■ ,XK~l) 



Oi 





(77) 
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where Q{xo, ■ ■ ■ ,xx/2-i) is a maximal rate square GLCOD 
of size N/2 [11], [19], x^ = Re{xi} + jIm{x(^,+K/2)K} 
where (a)x denotes a (mod K). 

Proof: The proof is by direct verification. As the maximal 
rate of square GLCOD of size N/2 is [1 1], [19] the rate 

of 5 in (|77) is 2^ 



and hence S is maximal rate. Next 
we show that 5 is a RFSDD. Consider 



S^S 











by construction, the sum of weight matrices of x^j , x^q for 
any symbol Xk is In and ( I44t -(I45> are satisfied as is a 
GLCOD. Therefore is a RFSDD. ■ 
Other square RFSDDs can be constructed from (I77> by 
applying some of the following 

permuting rows and/or columns of ( I77> . 

permuting the real symbols {xki,XkQ}^SQ, 

multiplying a symbol by -1 or ±j 

conjugating a symbol in J77t . 
Following [11, Theorem 2] we have 

Theorem 24: All square RFSDDs can be constructed from 
RFSDD S of (I77> by possibly deleting rows from a matrix of 
the form 

S' = USV (78) 

where U, V are unitary matrices, up to permutations and 
possibly sign change in the set of real and imaginary parts 
of the symbols. 

Proof: This follows from the observation after i69\ that 
the pair of sets {C'j.^}^S^, i — 1,2 constitute a Hurwitz family 
and Theorem 2 of [11] which applies to Hurwitz families. ■ 
It follows that the CIODs presented in Example l4.2l are unique 
up to multiplication by unitary matrices. Moreover, observe 
that the square RFSDDs of Theorem |23l can be thought of as 
designs combining co-ordinate interleaving and GLCODs. We 
therefore, include such RFSDDs in the class of co-ordinate 
interleaved orthogonal designs (CIODs), studied in detail in 
the next section. 

VI. Co-ordinate Interleaved Orthogonal Designs 

In the Section II VI we characterized SDDs in terms of the 
weight matrices. Among these we characterized a class of full- 
rank SDD called FSDD and classified it into UFSDD and 
RFSDD. In the previous section we derived and constructed 
maximal rate FSDDs. However, we have not been able to 
derive the coding gain of the either the class SDD or FSDD 
in general; the coding gain of GLCODs is well-known. This 
section is devoted to an interesting class of RFSDD C FSDD 
called Co-ordinate Interleaved Orthogonal Designs (CIODs) 
for which we will not only be able to derive the coding gain 
but also the Maximum Mutual Information. 

We first give an intuitive construction of the CIOD for two 
transmit antennas and then formally define the class of Co- 
ordinate Interleaved Orthogonal Designs (CIODs) comprising 
of only symmetric designs and its generalization. Generalized 
CIOD (GCIOD) which includes both symmetric and non- 
symmetric (as special cases) designs in Sub-section IVI-AI 
Also, we show that rate-one GCIODs exist for 2, 3 and 4 



transmit antennas and for all other antenna configurations the 
rate is strictly less than 1. A construction of GCIOD is then 
presented which results in rate 6/7 designs for 5 and 6 transmit 
antennas, rate 4/5 designs for 7 and 8 transmit antennas and 
rate ^^^ti GCIOD for 2m - 3, 2to - 2 > 8 coiTespo nd- 
ing to whether N is odd or even. In Subsection IVI-A.2I the 
signal set expansion associated with the use of STBC from 
any co-ordinate interleaving when the uninterleaved complex 
variables take values from a signal set is highlighted and the 
notion of co-ordinate product distance (CPD) is discussed. 
The coding gain aspects of the STBC from CIODs constitute 
Subsection IVI-BI and we show that, for lattice constellations, 
GCIODs have higher coding gain as compared to GLCODs. 
Simulation results are presented in Subsection IVI-CI The 
Maximum Mutual Information (MMI) of GCIODs is discussed 
in Subsection I VI-DI and is compared with that of GLCODs to 
show that, except for = 2, CIODs have higher MMI. In 
a nutshell this section shows that, except for A^ = 2 (the 
Alamouti code), CIODs are better than GLCODs in terms 
of rate, coding gain, MMI and BER. 

A. Co-ordinate Interleaved Orthogonal Designs 

We begin from an intuitive construction of the CIOD for two 
transmit antennas before giving a formal definition (Definition 
0. Consider the Alamouti code 



Xo Xi 



-X. 



X, 



When the number of receive antennas M — 1, observe that the 
diversity gain in the Alamouti code is due to the fact that each 
symbol sees two different channels ho and hi and the low ML 
decoding complexity is due to the use of the orthogonality of 
columns of signal transmission matrix, by the receiver, over 
two symbol periods to form an estimate of each symbol. 

Alternately, diversity gain may still be achieved by trans- 
mitting quadrature components of each symbol separately on 
different antennas. More explicitly, consider that the in-phase 
component, xqj, of a symbol, xq = a;o/+ja;oQ, is transmitted 
on antenna zero and in the next symbol interval the quadrature 
component, xqq, is transmitted from antenna one as shown in 
Table m 

It is apparent that this procedure is similar to that of co- 
ordinate interleaving (see Remark^] for references) and that 
the symbol has diversity two if the difference of the in-phase 
and quadrature components is not-zero, but the rate is half. 
This loss of rate can be compensated by choosing two symbols 
and exchanging their quadrature components so that one co- 
ordinate of each symbol is transmitted on one of the antennas 
as shown in Table HH 

As only one antenna is used at a time for transmission, the 
only operation required at the receiver to decouple the symbols 
is to exchange the quadrature components of the received 
signals for two symbol periods after phase compensation. 

The CIOD for four antennas is linked to the CIOD for two 
antennas in a simple manner The CIOD for two antennas uses 
complex symbols and uses antenna cycling between antennas 
and 1 . For four antennas consider antennas and 1 as one set 
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and antennas 2 and 3 as another set. Using two antennas and 
complex symbols, we can transmit a quaternion symbol (four 
co-ordinates) rather than a complex symbol (two co-ordinates). 
After interleaving the co-ordinates of the quaternion symbol 
we cycle between the first and second set of antennas. 

That the decoding is single-symbol decoding with the in- 
phase and quadrature-phase components having got affected 
by noise components of different variances for any GCIOD 
is shown in Subsection IVI-A.ll In the same subsection the 
fuU-rankness of GCIOD is also proved. If we combine, the 
Alamouti scheme with co-ordinate interleaving we have the 
scheme for 4 transmit antennas of Example I4.2l and whose 
receiver structure is explained in detail in Example 16.21 Now, 
a formal definition of GCIODs follows: 

Definition 7 (GCIOD): A Generalized Co-ordinate Inter- 
leaved Orthogonal Design (GCIOD) of size A^i x N2 in 
variables Xi,i = 0, ■ ■ ■ ,K — 1 (where K is even) is a L x 
matrix S{xq, - ■ ■ , xk-i), such that 



9i(io, 











(79) 



where Oi{xo, ■ ■ ■ , Xk/2-i) and 02ixK/2, ■ ■ ■ , xk-i) are 
GLCODs of size Li x Ni and L2 x respectively, with 
rates K/2Li, K/2L2 respectively, where iVi + iV2 = N, 
Li + L2 = L, Xi = Re{xi} + jIm{x(^i+K/2)K} and {a)K 
denotes a (mod K). If Qi = Q2 then we call this design a 
Co-ordinate interleaved orthogonal design(CIOD)^. 
Naturally, the theory of CIODs is simpler as compared to 
that of GCIOD. Note that when 61 = 62 and = L we 
have the construction of square RFSDDs given in Theorem 
I23I Examples of square CIOD for iV = 2,4 were presented 
in Example 14.21 

Example 6.1: An example of GCIOD, where Oi 7^ G2 is 
S{xo,-- ■ ,X3) 



XQI +iX2Q 
-Xii +jX3Q 






Xll 

xoi 



'3X3Q 
'jX2Q 



X21 ^ 



3X0Q 
\-jxiQ 



(80) 



where 81 is the rate-one Alamouti code and 82 is the trivial, 
rate-one, GLCOD for iV = 1 given by 



82 - 



Xq 

-xt 



Observe that S is non-square and rate-one. This code can also 
be thought of as being obtained by dropping the last column 
of the CIOD in J41> . Finally, observe that ( I80t is not unique 
and we have different designs as we take 



82 



etc. for the second GLCOD. 

'These designs were named as Co-ordinate interleaved orthogonal design 
(CIOD) in [47], [48] since two different columns are indeed orthogonal. 
However, the standard dot product of different columns may be different 
whereas in conventional GLCODs apart from orthogonality for two different 
columns, all the columns will have the same dot product. 



Xq 




Xo 


Xl_ 


1 


_-xi_ 



1 ) Coding and Decoding for STBCs from GCIODs: First, 
we show that every GCIOD is a RFSDD and hence is SD and 
achieves full diversity if the indeterminates take values from 
a signal set with non-zero CPD. 

Theorem 25: Every GCIOD is an RFSDD. 

Proof: Let be a GCIOD defined in We have 

8l^8i 
8^82 

a-klNt 
bklN2\ 



S^S = 



(81) 



(82) 



where Ofc = {j2k=o 



^2 

kl 



(k+K/2)KO I ^"'^ 

Observe that there are no 



\Y^k=K/2 Hi + ^lk+K/2)KQ 

terms of the form XkiXkQ,XkixiQ etc. in S'^S, and therefore 
is a SDD (this is clear from i22\ ). Moreover, by construe 



tion, the sum of weight matrices of x^j and x^q for any 
symbol x^ is /jv and hence 5* is a FSDD. Furthermore, for 
any given fc,0 < k < K — 1 the weight matrices of both 
x\j, x\q are not full-rank and therefore, by Definition|6j S is 
a RFSDD. ■ 
The transmission scheme for a GCIOD, S{xq, • • • , xk-i) 
of size N, is as follows: let Kb bits arrive at the encoder in 
a given time slot. The encoder selects K complex symbols. 
Si, i = 0, • • • , X — 1 from a complex constellation A of size 
1^1 = 2^ . Then setting Xi — Si^i ~ Q, - ■ ■ , if — 1, the encoder 
populates the transmission matrix with the complex symbols 
for the corresponding number of transmit antennas. The cor- 
responding transmission matrix is given by ^(so, • • • , sk-i)- 
The received signal matrix (|5} is given by. 



V = 5H + W. 



(83) 



Now as every GCIOD is a RFSDD (Theorem|25}, it is SD and 
the receiver uses (12 1> to form an estimate of each Si resulting 
in the ML rule for each Si, i = 0, • • • , A' — 1, given by 



min Mi{si) 



min II V - {A2iSii + yl2i+iSiQ)H|| 



(84) 



Remark 26: Note that forming the ML metric for each vari- 
able in ( I84L implicitly involves co-ordinate de-interleaving, in 
the same way as the coding involves co-ordinate interleaving. 
Also notice that the components sn and s^g (i.e., the weight 
matrices that are not full-rank) have been weighted differently 
- something that does not happen for GLCODs. We elaborate 
these aspects of decoding GCIODs by considering the decod- 
ing of rate-one, CIOD for 4 in detail. 

Example 6.2 (Coding and Decoding for CIO D for N = A): 
Consider the CIOD for = 4 given in i41\ . If the signals 
sq, si, S2, S3 € A are to be communicated, their interleaved 
version as given in Definition are transmitted. The signal 
transmission matrix, S, 

S0I+3S2Q S1I+]S3Q 



S = 



SQ 31 

-Sl/+jS3Q So/-jS2Q 

S2I+jSoQ Ssi+jSiQ 



-S3I+JS1Q S2I - JSOQ 



(85) 
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is obtained by replacing Xi in the CIOD by s; where each 
Si, i = 0,1,2,3 takes values from a signal set A with 2^ 
points. 

The received signals at the different time slots, Vjt, t — 



0,1,2,3 and j 
are given by 



0, 1, • • • , M — 1 for the M receive antennas 



VjO 
Vj3 



= hojSo + hijSi + Ujo; 

= /l2jS2 + ^3iS3 + nj2; 
= —h2jsl + /l3j,S2 + nj3 



where nji ,i — 0, 1, 2, 3 and j = 0, ■ 



(86) 

, Af — 1 are complex 



independent Gaussian random variables. 



Let V,- 



VjO, Vji, Vj2, 



= [fijo, nji, nj2, 



H, 



and 



S 



[So, 



si, S2, sa]^, 



hoj 


























h2j 










h* 


— h* 



where j = 0, 1, 
written as 



.M- 
V., 



1. Using this notation, ( I86> can be 
Wj. (87) 



Let 



VjO, Vjl, Vj2, Vj3 



Then, we have 



-H«W,. 



{\h 



2jl 



S 

(88) 



Rearranging the in-phase and quadrature-phase components 
of Vji's, (which corresponds to deinterleaving) define, for i — 
0,1, 



Vi 



M-l 
M-l 



ivji+2,Q = asij +ibsi^Q + uoi (89) 



bs^+2.i +iasi+2,Q + uii (90) 



where a = Ejlo^ilNP 



|/i3jp} and uoi,uii are complex Gaussian random variables. 
Let Wj = [ujQ fiji nj2 fij^]^ ~ H^Wj. Then uqi = 

j=o "j^J + i^ji+2.Q and = Z^^^q "^ji+^J + J"-r*,Q 
where i — 0, 1. Note that uoo and uoi have the same 
variance and similarly uio and uu. The variance of the in- 
phase component of uoo is a and that of the quadrature-phase 
component is b. The in-phase component of uio has the same 
variance as that of the quadrature-phase component of uqq and 
vice versa. The ML decision rule for such a situation, derived 
in a general setting is: Consider the received signal r, given 
by 

r = cisj +jc2SQ +71 (91) 

where ci,C2 are real constants and s/, sq are in-phase and 
quadrature-phase components of transmitted signal s. The 



ML decision rule when the in-phase, n/, and quadrature- 
phase component, rig, of the Gaussian noise, n have different 
variances cia^ and C2cr^ is derived by considering the pdf of 
n, given by 



Pn{n) 



1 



2-Ka'^^ciC2 



(92) 



The ML rule is: decide in favor of Si, if and only if 

Pn{r/si) > pn{r/sk), y i ^ k. (93) 

Substituting from ( 19 U and (|92} into ( I93> and simplifying we 
have 

C2k/-aSi,/P + Ci\rQ - bsi^ql'^ < C2\ri - askjl'^ 

+ci\rQ - bsk,Q\^ V^^fc. (94) 

We use this by substituting ci — a and C2 — b, to obtain j95t 
and ci — b and C2 = a, to obtain ( I96> . For Uj, j = 0,1, 
choose signal Si G A iff 

+ a|wj,Q - &Sfc,Qp, yi^k (95) 

and for Uj , j = 2, 3, choose signal Si iff 

a\vjj - bsijl^ + b\vj^Q - asi^gl'^ < a\vjj - bskjl'^ 

+b\vj.,Q - ask,Qf, Vi^k. (96) 

From the above two equations it is clear that decoupling of the 
variables is achieved by involving the de-interleaving operation 
at the receiver in ( I89> and ( I90l l. Remember that the entire 
decoding operation given in this example is equivalent to using 
j84l l. We have given this example only to bring out the de- 
interleaving operation involved in the decoding of GCIODs. 

Next we show that rate-one, GCIODs (and hence CIODs) 
exist for TV = 2, 3, 4 only. 

Theorem 27: A rate-one, GCIOD exists iff TV = 2, 3, 4. 
Proof: First observe from il9i that the GCIOD is rate- 
one iff the GLCODs 8i, 02 are rate-one. Following, Theorem 
|6] we have that a rate-one non-trivial GLCOD exist iff N — 2. 
Including the trivial GLCOD for = 1, we have that rate-one 
GCIOD exists iff iV = 1 + 1, 1 + 2, 2 + 2, i.e. iV = 2, 3, 4. ■ 
Next we construct GCIODs of rate greater than 1/2 for TV > 4. 
Using the rate 3/4 GLCOD i.e. by substituting 6i = 62 by the 
rate 3/4 GLCOD in (|79}, we have rate 3/4 CIOD for 8 transmit 
antennas which is given in ( I97l i. Deleting one, two and three 
columns from S we have rate 3/4 GCIODs for N — 7, 6, 5 
respectively. Observe that by dropping columns of a CIOD we 
get GCIODs and not CIODs. But the GCIODs for = 5, 6, 7 
are not maximal rate designs that can be constructed from the 
Definitional using known GLCODs. 

Towards constructing higher rate GCIODs for A^ = 5,6, 7, 
observe that the number of indeterminates of GLCODs 6i, 02 
in Definition are equal. This is necessary for full-diversity 
so that the in-phase or the quadrature component of each 
indeterminate, each seeing a different channel, together see all 
the channels. The construction of such GLCODs for A^i ^ N2, 
in general, is not immediate. One way is to set some of 
the indeterminates in the GLCOD with higher number of 
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Six), 



S4{X3I +ixoQ,X4i +ixxQ,X5i +}X2q) 



(97) 



indeterminates to zero, but this results in loss of rate. We next 
give the construction of such GLCODs which does not result 
in loss of rate. 

Construction 6.1: Let 8i be a GLCOD of size Li x 
Ni, rate ri — Ki/Li in Ki indeterminates xq,--- ,xki-i 
and similarly let 82 be a GLCOD of size L2 x N2, rate 
r2 — K2/L2 in K2 indeterminates yo, - ■ ■ ,yK2-i- Let K = 
lcm{Ki, K2), rii — K/Ki and n2 = K/K2- Construct 



81 



and 



82 - 



8i(a;o,a;i,- • • ,XKt~i) 

Ql{XKt , XKi + l, ■ ■ ■ , X2Kt-l) 
Ql{x2Ki,X2Kt + l, - ■ ■ ,XzKi-l) 

2:(ni-l)ifi + l, • • • 1 XniKi-l) 



62(2/0,2/1, • • • ,2/i^2-l) 
©2(2/^^2, 2//f2+i, •• • ,2/2^2-1) 

©2(2/2^^2,2/2^2 + 1,- • • ,2/3if2-l) 
02(2/(n2-l)K2' 2^(n2-l)-fs:2 + l' ■ ■ ■ ,2/n2if2-l) 



(98) 



(99) 

Then 81 of size niLi x A^i is a GLCOD in indeterminates 
xo, xi, • • • , xk-1 and 82 of size rt2-^2 x Af2 is a GLCOD in 
indeterminates 2/0,2/1, • • • ,2/if-i- Substituting these GLCODs 
in (|79j we have a GCIOD of rate 



n = 



2K 



niLi + 712X2 
2\cm{Ki,K2) 



niLi + 712^2 

2\cm{Ki,K2) 

lcm{Ki,K2){Li/Ki + L2/K2) 
= Hin,r2) (100) 

where H{ri,r2) is the Harmonic mean of ri,r2 with N = 
A^i + N2 and delay, L = niLi + 7i2i2- 
We illustrate the Construction 16. II bv constructing a rate 6/7 
GCIOD for six transmit antennas in the following example. 
Example 6.3: Let 



81 



.To Xi 

— Xi Xq 



be the Alamouti code. Then Li = Ni = Ki = 2. Similarly 
let 



Then L2 



82 



N2 



Xq 



Xl X2 
Xl X2 

Xq —Xi 



Xq 



A, K2 = 3 and the rate is 3/4. K 



lcm{Ki,K2) = 6,ni = K/Ki = 3 and = K/K2 = 2. 



81 = 



8i(a;o,a;i) 

0i(a;2,a:3) 
8i(a;4,a:5) 



Xq 


Xl 


-xl 


Xq 


X2 


X3 




X2 


X4 


X5 






-4 





(101) 



Similarly, 



82 = 



Xq 

-xl 

— X2 



X3 







The GCIOD for N ^ 
rate of the GCIOD in 



Ni - 



Xl 

Xq 



—X2 

X4 



-N2 = 
is = 

14 



X2 



Xq 

Xz 







X2 
-Xl 
Xq 



x^ 

— Xi 

X3 



(102) 



6 is given in ( I103t . The 
f 0.8571 > 3/4. This 
increased rate comes at the cost of additional delay. While the 
rate 3/4 CIOD for = 6 has a delay of 8 symbol durations, 
the rate 6/7 GCIOD has a delay of 14 symbol durations. In 
other words, the rate 3 /4 scheme is delay-efficient, while the 
rate 6/7 scheme is rate-efficient**. Deleting one of the columns 
we have a rate 6/7 design for 5 transmit antennas. 
Similarly, taking 81 to be the Alamouti code and 82 to be the 
rate 2/3 design of [17] in Construction 16. II we have a CIOD 
for N = 7 whose rate is given by 



7^ = 



3/2 + 1 



= -.0.8. 



n + 2 
= n/2 



We have the following theorem: 

Tlieorem 28: The maximal rate of GCIOD for A^ 
antennas, TZ is lower bounded as 7?. > ^i"'"!"^''' where r 
if n is even or to = {n + l)/2 if n is odd. 

Proof: We need to prove that a GCIOD of rate TZ > 
^3m+i'' where m — n/2 if n is even or m = [n + l)/2 if n 
is odd exists. 

Consider Construction 16.11 For a given A^, Let 81 be the 
Alamouti code. Then Li = A^i = ATi = 2 and Ar2 = A^ - 2. 
Let 82 be the GLPCOD for n — N — 2 transmit antennas with 
rate r2 — where m = n/2 if n is even or m — {n + l)/2 
if n is odd [17]. The corresponding rate of the GCIOD is given 
by 

2 _2(to+1) 



7^ = 



2m 
m+1 



+ 1 



3to + 1 



Significantly, tliere exist CIOD and GCIOD of rate greater 
tliat 3/4 and less than 1, while no such GLCOD is known 
to exist. Moreover for different choice of 81 and 82 we have 
GCIODs of different rates. For example: 

^Observe that we are not in a position to comment on the optimaUty of 
both the delay and the rate. 
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Xll 


+ j^rQ 




















-Xil 


+ ix7Q 


^01 






















X2I - 






+ i^9Q 




















-X3I 


+ }X9Q 


























-j^lOQ 


























+ 


X41 


- j^lOQ 
































a; 71 






+ 'iX2Q 

















-xri +}XiQ 


a;6l 











'r3X2Q 












-XgJ +iX2Q 







X6I 


- i^OQ 


-X-JI 

















-^ai 








Xdl - 
















^101 




^111 


+ 3^bQ 

















-OSlQI +jX4Q 


3^91 























-Xiil +}X5Q 







X9I 




-^101 



















+ jx5Q 




- }^4Q 


^91 - 


'r'iXiQ _ 



(103) 



Example 6.4: For a given N, Let Qi be the Alamouti code. 
Then Li= Ni^ Ki^2 and iVa = Let 62 be the rate 

1/2 GLPCOD for N — 2 transmit antennas (either using the 
construction of [9] or [15]). Then r2 — 1/2. The corresponding 
rate of the GCIOD is given by 



In Table mil we present the rate comparison between GLCODs 
and CIODs-both rate-efficient and delay efficient; and in Table 
II VI we present the delay comparison. 

Observe that both in terms of delay and rate GCIODs are 
superior to GLCOD. 

2) GCIODs vs. GLCODs: In this subsection we summarize 
the differences between the GCIODs and GLCODs with 
respect to different aspects including signal set expansion, 
orthogonality and peak to average power ratio (PAPR). Other 
aspects like coding gain, performance comparison using simu- 
lation results and maximum mutual information are presented 
in subsequent sections. 

As observed earlier, a STBC is obtained from the GCIOD by 
replacing Xi by Si and allowing each Si, i = 0,1, ■ ■ ■ , K — 1, 
to take values from a signal set A. For notational simplicity we 
will use only S for S{xo, ■ ■ ■ , xk-i) dropping the arguments, 
whenever they are clear from the context. 

The following list highlights and compares the salient 
features of GCIODs and GLCODs: 

. Both GCIOD and GLCOD are FSDD and hence STBCs 
from these designs are SD. 

. GCIOD is a RFSDD and hence STBCs from GCIODs 
achieve full-diversity iff CPD of A is not equal to zero. 
In contrast STBCs from GLCODs achieve full-diversity 
for all A. 

. Signal Set Expansion: For STBCs from GCIODs, it 
is important to note that when the variables Xi, i — 
0,1, ■ ■ ■ , K ~ 1, take values from a complex signal set 
A the transmission matrix have entries which are co- 
ordinate interleaved versions of the variables and hence 
the actual signal points transmitted are not from A but 
from an expanded version of A which we denote by A. 
Figure |2ja) shows A when A = {!,— l,j,— j} which 
is shown in Figure |2jc). Notice that A has 8 signal 
points whereas A has 4. Figure |3b) shows A' where 
A' is the four point signal set obtained by rotating 
A by 13.2825 degrees counter clockwise i.e.. A' ~ 
{ej^-ej^jej^-jeJ'^} where 6 = 13.2825 degrees as 



shown in Figure |3d). Notice that now the expanded 
signal set has 16 signal points (The value 9 = 13.2825 
has been chosen so as to maximize the parameter called 
Co-ordinate Product Distance of the signal set which is 
related to diversity and coding gain of the STBCs from 
GCIODs, discussed in detail in Section lVI-Br It is easily 
seen that \A!\ < \A\'^. 

Now for GLCOD, there is an expansion of signal set, 
but \A'\ < 2\A\. For example consider the Alamouti 
scheme, for the first time interval the symbols are from 
the signal set A and for the next time interval symbols 
are from A*, the conjugate of symbols of A. But for 
constellations derived from the square lattice \A'\ << 
2\A\ and in particular for square QAM \A'\ = \A\. So 
the transmission is from a larger signal set for GCIODs 
as compared to GLCODs. 

• Another important aspect to notice is that for GCIODs, 
during the first L/2 time intervals iVi < iV of the N 
antennas transmit and the remaining iV2 = N — Ni 
antennas transmit nothing and vice versa. So, on an 
average half of transmit antennas are idle. 

• For GCIODs, S, is not an scaled orthonormal matrix but 
is an orthogonal matrix while for square GLCODs, S, is 
scaled orthonormal. For example when S is the CIOD 
given by ( I85t for = 4 transmit antennas, 

. GCIODs out perform GLCODs for TV > 2 both in terms 
of rate and delay as shown in Tables |lll| and IIVI 

• Due to the fact that at least half of the entries of GCIOD 
are zero, the peak-to-average power ratio for any one 
antenna is high compared to those STBCs obtained from 
GLCODs. This can be taken care of by "power uni- 
formization" techniques as discussed in [11] for GLCODs 
with some zero entries. 



B. Coding Gain and Co-ordinate Product Distance (CPD) 

In this section we derive the conditions under which the 
coding gain of the STBCs from GCIODs is maximized. 
Recollect from Section HVl that since GCIOD and CIOD are 
RFSDDs, they achieve full-diversity iff CPD of A is non- 
zero. Here, in Subsection IVI-B.ll we show that the coding 
gain defined in is equal to a quantity, which we call, 
the Generalized CPD (GCPD) which is a generalization of 
CPD. In Subsection IVI-B.21 we maximize the CPD for lattice 
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|ioP + |xiP 








|ioP + |xiP 








\x2? + \xi\^ 









|i2p + |i3| = 



(104) 



constellations by rotating the constellation''. Similar results are 
also obtained for the GCPD for some particular cases. We 
then compare the coding gains of STBCs from both GCIODs 
and GLCODs in Subsection IVI-B.5I and show that, except 
for N = 2, GCIODs have higher coding gain as compared 
to GLCODs for lattice constellations at the same spectral 
efficiency in bits/sec/Hz. 

1 ) Coding Gain of GCIODs: Without loss of generality, 
we assume that the GLCODs Qi, 82 of Definition are such 
that their weight matrices are unitary. Towards obtaining an 
expression for the coding gain of CIODs, we first introduce 

Definition 8 (Generalized Co-ordinate Product Distance): 
For arbitrary positive integers A^i and the Generalized 
Co-ordinate Product Distance (GCPD) between any two 
signal points u = uj + juq and v ^ vj + jwg, u 7^ w of 
the signal set A is defined in (I1Q5> and the minimum of this 
value among all possible pairs of distinct signal points of the 
signal set A is defined as the GCPD of the signal set and will 
be denoted by GCPi^Ar^^jVa (-4) or simply by GCPDn^^n^ 
when the signal set under consideration is clear from the 
context. 

Remark 29: Observe that 

1) When Ni= N2, the GCPD reduces to the CPD defined 
in Definition |4] and is independent of both A^i and N2 . 

2) GCPDmi,N2{'^tv) = GCPDn2,Ni{u,v) for any two 
signal points u and v and hence GCPDmi,N2{^) = 

GGPDm^M^)- 
We have. 

Theorem 30: The coding gain of a full-rank GCIOD with 
the variables taking values from a signal set, is equal to the 
GCPI? Ari,jV2 of that signal set. 

Proof: For a GCIOD in Definition we have. 



S^S = 



aKlNi 

bKlN2 

|ioP H \- \xk/2-i?, 



(106) 



\Xk/2\ 



\XK-l\ 



where Xi = Re{xi} + j/TO{a;(i+K/2)A-} ™d where (a)^ 
denotes a (mod K). Consider the codeword difference matrix 
B{S, S') = S - S' which is of full-rank for two distinct 
codeword matrices S, S'. We have 



B^(S,S')S(S,S') = 



V^k/at^ 



z.' |2 



VflA' = |Xo-.TorH ^ I^K/2-1 - 5^/2-1 1 ): 



\xk/2-S:k/2\ ^ ^ |iA'-l 



'The optimal rotation for 2-D QAM signal sets is derived in [62] using 
Number theory and Lattice theory. Our proof is simple and does not require 
mathematical tools from Number theory or Lattice theory. 



where at least one Xk differs from x'j,, k = 0, • • • ,K— 1. 
Clearly, the terms (|io — a;oP + ' • • + 1^x72-1 —ix/2-iP) ^"'^ 
(|iif/2~i/f/2p + ' ■ • + |a5A'-i — are both minimum iff 
Xk differs from x'/, for only one k. Therefore assume, without 
loss of generality, that the codeword matrices S and S' are such 
that they differ by only one variable, say xq taking different 
values from the signal set A. Then, for this case. 



Ai 



= I To/ 



det{B^(S,S')-B(S,S')} 

X'0I\^^\XQQ- 



1/N 



2N2 



Similarly, when S and S' are such that they differ by only in 
xk/2 then 

A2 = det{B^(S,S')B(S,S')}'^^ 

2Ni 

= \xk/2I-x'k/2i\'''^\xk/2Q~x'k/2q\^'^ 

and the coding gain is given by minjjo^xjc/^e-A {Ai, A2} = 

GGPDn,,N2- ■ 

An important implication of the above result is. 

Corollary 31: The coding gain of a full-rank STBC from 
a CIOD with the variables taking values from a signal set, is 
equal to the CPD of that signal set. 

Remark 32: Observe that the CPD is independent of the 
parameters A^i , N2 and is dependent only on the elements of 
the signal set. Therefore the coding gain of STBC from CIOD 
is independent of the CIOD. In contrast, for GCIOD the coding 
gain is a function of A^i, N2. 

The full-rank condition of RFSDD i.e. GPD ^ can be 
restated for GCIOD as 

Theorem 33: The STBC from GCIOD with variables tak- 
ing values from a signal set achieves full-diversity iff the 
GCPD]y-^ N2 of that signal set is non-zero. 
It is important to note that the GGPDni,N2 is non-zero iff 
the CPD is non-zero and consequently, this is not at all a 
restrictive condition, since given any signal set A, one 
can always get the above condition satisfied by rotating it. 
In fact, there are infinitely many angles of rotations that 
will satisfy the required condition and only finitely many 
which will not. Moreover, appropriate rotation leads to 
more coding gain also. From this observation it follows that 
signal constellations with GPD = and hence GGPD ~ 
like regular M — ary QAM, symmetric M — ary PSK will 
not achieve full-diversity. But the situation gets salvaged by 
simply rotating the signal set to get this condition satisfied 
as also indicated in [42], [43], [60]. This result is similar to 
the ones on co-ordinate interleaved schemes like co-ordinate 
interleaved trellis coded modulation [42], [43] and bit and 
co-ordinate interleaved coded modulation [40]-[45], [55] for 
single antenna transmit systems. 

2) Maximizing CPD and GCPD for Integer Lattice con- 
stellations: In this subsection we derive the optimal angle of 
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GCPDn,,n2{u,^) 



2N2 



2N2 



(105) 



rotation for QAM constellation so that the CPD and hence 
the coding gain of CIOD is maximized. We then generalize 
the derivation so as to present a method to maximize the 
GCPDn, 

3) Maximizing CPD: In the previous section we showed 
that the coding gain of CIOD is equal to the CPD and 
that constellations with non-zero CPD can be obtained by 
rotating the constellations with zero CPD. Here we obtain the 
optimal angle of rotation for lattice constellations analytically. 
It is noteworthy that the optimal performance of co-ordinate 
interleaved TCM for the 2-D QAM constellations considered 
[42], [43], using simulation results was observed at 32°; 
analytically, the optimal angle of rotation derived herein is 
9 = tan(2)/2 = 31.7175° for 2-D QAM constellations. The 
error is probably due to the incremental angle being greater 
than or equal to 0.5. We first derive the result for square QAM. 

Theorem 34: Consider a square QAM constellation A, 
with signal points from the square lattice {2k — 1 — Q)d + 
j(2Z — 1 — Q)d where k,l e [1,Q] and d is chosen so that 
the average energy of the QAM constellation is 1. Let 9 be 
the angle of rotation. The maximum CPD of A is obtained 
at 9opt = ^'''^^^"^-^^ — 31.7175° and is given by 

CPDopt - (108) 

Proof: The proof is in three steps. First we derive 
the optimum value of 9 for 4-QAM, denoted as 9opt (the 
corresponding CPD is denoted as CPDopt)- Second, we 
show that at 9opt, CPDopt is in-fact the CPD for all other 
(square) QAM. Finally, we show that for any other value of 
6* e [0,7r/2], CPD < CPDopt completing the proof. 
Step 1: Any point P(x, y) e rotated by an angle 9 e 
[0, 90°] can be written as 

cos 9 sin 9 
— sin 9 cos 9 







VR 





X 

y 



(109) 



Let Pi{xi,yi), P2{x2,y2) be two distinct points in A such 
that Ax — xi ~ X2, Ay = j/i — y2- Observe that Ax, Ay = 
0, ±2d, ■■■ , ±2(Q - l)d. We may write Ax = ±2md, Ay = 
±2nd, m, n G [0,Q ~ 1] but both Ax. Ay cannot be zero 
simultaneously, as Pi , P2 are distinct points in A- Since, 
rotation is a linear operation. 



AXr 

Ayr 



R 



Ax 
Ay 



(110) 



where Ax^ = Xi^ - X2R,Ayr = yi^ - y2R- The 
CPD between points Pi and P2 after rotation, denoted by 
CPD{Pir, P2r), is then given by 



CPD{Pir,P2r) = |Ax,||Ay,| = 



AxAycos(26') 



(Ax)2 - (Ay)2 



sin(26l) 



(111) 



For 4-QAM, possible values of CPD{Pir, P2r) are 
CPDi{Pir,P2r) = 2d^\sm{29)l 

CPD2{Plr,P2r) = 4^^ | COs(20) | . 



(112) 



Fig. |3] shows the plots of both CPDi and CPD2- As sine 
is an increasing function and cosine a decreasing function 
of 9 in the first quadrant, equating CPDi,CPD2 gives the 
optimal angle of rotation, 9opt- Let CPD{9) be the CPD 
at angle 9 and CPDopt = ma.xg CPD{9)- It follows that 
9opt ^ '''■<=iM±^) ^ 31.7175°, 58.285° and CPDopt = 
2d^ sm{29opt) = M^cos{29opt) = ^- 

Step 2: Substituting the optimal values of 
sin(26'opt), cos(20opt) in il 1 1> we have for any two 
arbitrary points of a square QAM constellation, 

4d2 



CPD{Pir,P2r) 



V5 



\±nm + n —m | where n, m e Z 

(113) 



and both n, m are not simultaneously zero and Z is the set of 
integers. It suffice to show that 



± nm 



'1 > 1 Vn, m 



provided both n, m are not simultaneously zero. We consider 
the ± case separately. We have 



m 

"2" 

m 













' i) 









{1 + V5} 



Similarly, 



I — nm + ri^ 



The quadratic equation in n, 
roots 



2 

m 



± nm 



■{1 

{1- 



- V5} 



has 



-{±l±\/5}. 



Since n, ?Ti G Z, | ± nm - 



E Z and is equal to zero 



only if n = 0, ^{±l±\/5}. Necessarily, \±nm + n'^ — m'^\ > 
1 for n,m G Z and both n, m are not simultaneously zero. 
Therefore 9opt and CPDopt continue to be the optimum values 
of angle and the CPD for any square QAM. 
Step 3: Next we prove that for all other values of 61 e [0, 
CPD{9) < CPDopt- To this end, observe that for any value 
of 9 other than 9opt either CPDi or CPD2 is less than 
CPDopt (see the attached plot of CPDi,CPD2 in Fig. OJi. 
It follows that 



CPD{9) < CPD, 



opt 



with equality iff 6 — 6, 



opt ■ 
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Observe that Theoreml34lhas application in all schemes where 
the performance depends on the CPD such as those in [49], 
[44], [45], [42], [43], etc. and the references therein. 
Remark: The 4 QAM constellation in Fig. |2jc) is a rotated 
version (45°) of the QAM signal set considered in Theorem 

M 

Next we generalize Theorem |34] to all integer lattice con- 
stellations obtainable from a square lattice. We first find 
constellations that have the same CPD as the square QAM 
of which it is a subset. Towards that end we define. 

Definition 9 (NILC ): A Non-reducible integer lattice con- 
stellation (NILC) is a finite subset of the square lattice, 
{2k)d + j{2l)d where fc, Z e Z, such that there exists at least 
a pair of signal points pi = {2ki)d + j(2Zi)d and p2 — 
(2/c2)rf + j(2/2)rf such that either \ki - /C2I = 1, |/i - ^2! = 
or|;i-;2| = l,|fci-A:2|=0. 
We have. 

Corollary 35: The CPD of a non-reducible integer lattice 
constellation. A, rotated by an angle 9, is maximized at 9 = 

arctaii(2) 



2 = 31.7175° and is given by 

CPD, 



Proof: Since ^ is a subset of an appropriate square QAM 
constellation, we immediately have from Theorem 1341 

(115) 



CPDopt > 



V5" 



We only need to prove the equality condition. The CPD 
between any two points NILC at 9opt is given by ( II 13> 

4d2 



CPD{Pi,P2) = ^ \±nm + n^ -m^\ 
v5 



where n,m £ Z. 

(116) 

Since for NILC there exists at least a pair of signal points 
pi = {2ki)d + i{2li)d and p2 = {2k2)d + j{2l2)d such that 
either |fci-A:2| = l,\h-l2\ =0or|Zi-;2| = 1, |fci-A:2| = 0, 
we have CPL»(pi,P2) = ■ 
In addition to the NILCs, the lattice constellations that are a 
proper subset of the scaled rectangular lattices, (4fc)d + j(2?)d 
and {2k)d+i{Al)d where k,l e 'Z have CPD equal to 

All other integer lattice constellations have CPD > 

4) Maximizing the GCPD of the QPSK signal set: To derive 
the optimal angles of rotation for maximizing the GCPD we 
consider only QPSK, since the optimal angle is not the same 
for any square QAM, as is the case with CPD. 

Theorem 36: Consider a QPSK constellation A, with sig- 
nal points {2k ~ 3)d + j{2l - 3)d where k,l G [1,2] and 
d = l/\/2, rotated by an angle 9 so as to maximize the 
GCPDmj^,N2- The CCPDni,N2{^) is maximized at 9 opt ~ 
arctan(a:o) where is the positive root of the equation 



1 

1 

X 



2Ni 



(1 



xf''^ = l 



(117) 



where A^i > N2 and the corresponding GCPDni,N2{^) is 



1+^0 



Proof: Following the same notations as in Step 1 of 
Theorem |33 we have 

|Ax,.|^i|Ay^|^= = \2dmcos{9) + 2dnsui{9)f^ 

\~2dmsa\{9) + 2dncos(6')|^fl.l8) 

The possible values of GCPDi^j^^ j^^-j{Pi, P2) are 

2Ni 

GCPDi = Ad^ \sm{9) - cos{9)\T'f^+T^- 

|sin(6') + cos(6')|^^TTivi (119) 

2Ni 

GCPD2 = 4^2 |sin(6i) + cos{9)\T^Tn^ 

2N2 

|sin(6') -cos(6')|^^TTivi (120) 

2Ni 2N2 

GCPD3 = 4^2 |sin(6')|^^rTivi |cos(6»)|^^'rTivi (12I) 

2Ni 2N2 

GCPDi = 4^2 |cos(6l)|T^TTivi |sin(6')|''TTivi , (122) 

Now by symmetry it is sufficient to consider 9 e [0,7r/4). 
In this range sin(0) < cos{9) < 1 and accordingly, if 
iVi > N2 then GCPD3 < GCPDi and similarly GCPDi < 
GCPD2. Equating GCPDi, GCPD3 gives the optimal angle 
of rotation, 9 opt- We have 



(114) GCPDi ^GCPD3 



2JV2 



ism{9opt) - cos(6lopt))"i+"^ {sm{9opt) + cos(6'opt)) "1+"= 



2N2 



ism{9opt)) icos{9opt)) "1+"^ 



2N2 



{l-cot{9opt))"^+"^ (l + tan(e'opt))"i+"2 =1. 

Substituting ta,n{9opt) = x we have that x is the root of (^Hji 
The GCPDi and hence the GCPD at this value is 



GCPDi = 4d^ \sinieopt) - cos{9opt)\ 



JV1+JV2 



\sin{eopt)+cos{9opt)\"^ 



Ad 



Ad-" 



2Ni 2N2 

2 (Xo ~ l)«l+«2 [xq + l)"i + "2 



1 + xf, 



_ NI+N2 



1 + 2^0 



(123) 



Table fVl gives the optimal angle of rotation for various values 
of N — Ni + N2 along with the normalized GCPDM1.N2 
{GCPDisi^^N^)IAd^). Observe that for any given N the coding 
gain is large if iVi, iV2 are of the same size i.e., nearly equal. 
Also observe that the optimal angle of rotation lies in the range 
(26.656, 31.7175] and the corresponding normalized coding 
gain varies from (0.2,0.4472]. 

Note that the infimum corresponds to the limit where 
iVi = A^, A^2 = and the maximum corresponds to iVi — 
N2 = N/2. Unfortunately, the optimal angle varies with the 
constellation size, unlike CPD. In the next proposition we find 
upper and lower bounds on GCPDM1.N2 for rotated lattice 
constellations. 

Proposition 37: The GCPDmi,N2 for rotated NILC is 
bounded as 

2N2 

CPDT^TTN^ < GCPDn,,n2 < CPD, N2 > Ni 
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with equality iff A^i = A^2- 

Proof: Let p, q be two signal points such that 

GCPDn.^N, = GCPDN,,N,{p,q). 



(124) 



When A^i = N2 or Ax = Aj/ there is nothing to prove as the 
inequality is satisfied. 

Therefore let A^i 7^ N2 and Ax 7^ Ay. When the signal 
points are from the square lattice (2fc)d + j(2^)d where fc, I e 
Z and d is chosen so that the average energy of the QAM 
constellation is 1, rotated by an angle 6 then 

( 2JV2 
GCPD(^N^^N2){P,q) = min||Aa;^|"i+"2|A?/r|"i+"2, 

|Axr|™i+™2 |A?/r|"i+"2 \ 

2Ni 

= |mcos(6') + nsin(6')| "i+"2 

2N2 

|-TOsin(6') + 7icos(6')| "i+"fl?5) 

where rn, n G Z. For a NILC the GCPDni,N2 is upper 
bounded by the GGPDn^,N2 for QPSK and is given by (ITBl . 
Now the root of (II 17> . xo, is such that xo G (0.5, 1) and 
iV2 > N/2 and we immediately have 



■2N2 



4d^ 



(1 



Xo 



(1 



(126) 



completing GCPD 



Ni,N2 



< CPD. For the second part 
observe that, for N2 > Ni, \mcos{6) + nsm{d)\'^^ < 
\m cos{9) + n sm{9)\^'^ as |tocos(6') + risin(6')| < 1. Sub- 
stituting this in (I125> we have the lower bound. ■ 
In Proposition 1371 if we use 6 — arctan(2) for rotating the 
NILC then the GCPD is bounded as 



opt 



< GGPDn,m2 < CPDopu N2 > iVi,(127) 



2N2 



< gcpdn,,n2 < [-^ ], N2> mm) 



Remark 38: It is clear from Table |V] and the above in- 
equalities on GCPD that the value of GCPD decreases as the 
QAM constellation size increases and also as the difference 
between Ni,N2 increases. Therefore, while Construction 16.11 
gives high-rate designs, the coding gain decreases for QAM 
constellations. 

5) Coding gain of GCIOD vs that of GLCOD: In this 
subsection we compare the coding gains of GCIOD and 
GLCOD for the same number of transmit antennas and the 
same spectral efficiency in bits/sec/Hz-for same total transmit 
power For the sake of simplicity we assume that both GCIOD 
and GLCOD use square QAM constellations. 

6} The number of transmit antennas N=2: The total trans- 
mit power constraint is given by tr (S'^S') = L = 2. If the 
signal set has unit average energy then the Alamouti code 
transmitted is 



V2 



Xo 

-xt 



Xl 
Xn 



where the multiplication factor is for power normalization. For 
the same average transmit power the rate-one CIOD is 



xoi +3X1C 






Xll +jXQQ 



Therefore the coding gain of the Alamouti code for NILC is 
given by ^ and that of CIOD is given by Theorem|33as 
Therefore the coding gain of the CIOD for N=2 is inferior 
to the Alamouti code by a factor of ^ = 2% ~ 0.894, 
which corresponds to a coding gain of 0.4 dB for the Alamouti 
code'". 

7) The number of transmit antennas N=4: The average 
transmit power constraint is given by tr (S'^S') ~ L — A. If 
the signal set has unit average energy then the rate 3/4 COD 
code transmitted is 



1 

73 



Xo Xl 

-xl Xo 

-X*2 

. -X^2 



X2 



Xq 

xl 



■ 

X2 
~Xi 
XQ 



where the multiplication factor is for power normalization. For 
the same average transmit power, the rate 1 CIOD is given in 
(I129> . If the rate 3/4 code uses a 2" square QAM and the 
rate 1 CIOD uses a 2t square QAM, then they have same 
spectral efficiency in bits/sec/Hz, and the possible values of 
n for realizable square constellations is n = 8i,i G Z+. Let 
di , c?2 be the values of d so that the average energy of 2" 
square QAM and 2~ square QAM is 1. Therefore the coding 
gain of rate 3/4 COD for NILC is given by Kcod = ^ 
and that of CIOD is given by Theorem 1341 as Aqiod ~ 
Using the fact that for unit average energy M-QAM square 



constellations d 



M-l 



, we have 



A 



8 



COD = 



(28« - 1) 



and AcioD = 



12 



V5(26» 



1) 



where i G Z^ 



for a spectral efficiency of 6i bits/sec/Hz. For i = 1,2,3 we 
have AcoD = 0.0314, 1.2207e-004, 4.7684e-007 and Aciod 
= 0.0422, 6.5517e-004, 1.0236e-005 respectively, correspond- 
ing to a coding gain of 1.29, 7.29, 13.318 dB for the CIOD 
code. Observe that in contrast to the coding gain for N ^ 2 
which is independent of the spectral efficiency, the coding gain 
for = 4 appreciates with spectral efficiency. 

8) The number of transmit antennas N=8: The total trans- 
mit power constraint is given by tr (^S'^S^ = L — 8. If the 
signal set has unit average energy then the rate 1/2 COD code 
has a multiplication factor of 1/2 and for the same transmit 
power, the rate 3/4 CIOD has a multiplication factor of l/\/3. 
If the rate 1/2 COD code uses a 2" square QAM and the 
rate 3/4 CIOD uses a 2 "2" square QAM, then they have same 
spectral efficiency in bits/sec/Hz, and the possible values of 
n for realizable square constellations is n = 4i,i G Z+. Let 
di , c?2 be the values of d so that the average energy of 2" 
square QAM and 2 2 square QAM is 1 . Therefore the coding 
gain of rate 1/2 COD for NILC is given by Aqod 
and that of CIOD is given by Theorem 1341 as Aqiod 



4di 



3V5' 

'"in Section IVIII we revisit tliese codes for their use in rapid-fading 
channels where we show that this loss of coding gain vanishes and the CIOD 
for TV = 2 is SD while the Alamouti code is not. 
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1 

71 



-Xil +}X-iQ 






xii +^X3Q 









X2I +jxoQ 

-X3I + jXlQ 






X3I +jxiQ 
0=21 -j^OQ 



(129) 



Using the fact that for unit average energy M-QAM square 
constellations d = ,/^, we have 



A, 



6 



COD 



(24* - 1; 



and AcioD 



where i G Z"*" 



for a spectral efficiency of 3i bits/sec/Hz. For i ~ 1,2,3 
we have Acod= 0.4, 0.0235, 0.0015 and Aciod = 0.4737, 
0.0563, 0.007 respectively, corresponding to a coding gain of 
0.734, 3.789, 6.788 dB for the CIOD code. Observe that as in 
the case of = 4 the coding gain appreciates with spectral 
efficiency. 

Next we compare the coding gains of some GCIODs. 

9) The number of transmit antennas N=3: Both the GCIOD 
and GCOD for = 3 is obtained from the = 4 codes 
by dropping one of the columns, consequently the rates and 
the total transmit power constraint are same as for = 4. 
Accordingly, the rate 3/4 GCOD code uses a 2" square QAM 
and the rate-one GCIOD uses a 2~ square QAM where n — 
8i,ieZ+. The coding gain for the rate 3/4 GCOD for NILC is 
given by Aqcod — and that of GCIOD is lower bounded 

by Proposition 1371 as Agciod > ( 5^ ) ■ Using the fact that 
for unit average energy M-QAM square constellations d — 



0] shows the BER performance for these schemes. Observe 
that the scheme of this paper outperforms rate 1/2 COD by 
3.0 dB, rate 3/4 COD by 1.3 dB and STBC-CR by 1.2 dB 
at Pb = 10~^. A comparison of the coding gain. A, of these 
schemes is given in tabular form in Table IVll 

For CIOD, Aciod = 0.4478 while for STBC-CR 
^STBC-CR = 0.5 at i? = 2 bits/sec/Hz, but still CIOD out- 
performs STBC-CR because the coding gain is derived on 
the basis of an upper bound. If we take into consideration 
the kissing number i.e. the number of codewords at the given 
minimum coding gain, then we clearly see that though STBC- 
CR has higher coding gain, it has more than double the kissing 
number of CIOD. The results for rest of the schemes are in 
accordance with their coding gains; 

101og,„| . ^^^"^ 1=3.5 



A, 



rate 1/2 COD 



and 



10 log 



10 



A^ 



Ac/OD 
ate 3/4 COD 



1.3. 



M-l 



we have 



A, 



GOOD 



(28' - 1) 



and Agciod > 



12 



V5(26' _ 1) 



where 



for a spectral efficiency of 6i bits/sec/Hz. For i ~ 1,2,3 
we have Agcod = 0.0314, 1.22076-*, 4.76846"^ and 
^GCiOD > 0.0147, 5.69e-'\ 2.22e~'^ respectively. 

Observe that at high spectral rates, even the lower bound is 
larger than the coding gain of GCOD. In practice, however, 
the GCIOD performs better than GCOD at all spectral rates. 

C. Simulation Results 

In this section we present simulation results for 4-QAM and 
16-QAM modulation over a quasi-static fading channel. The 
fading is assumed to be constant over a fade length of 120 
symbol durations. 

First, we compare the CIOD for A^ = 4, with (i) the STBC 
(denoted by STBC-CR in Fig.Eland|5} of [62], (ii) rate 1/2, 
COD and (iii) rate 3/4 COD for four transmit antennas for 
the identical throughput of 2 bits/sec/Hz. For CIOD the 
transmitter chooses symbols from a QPSK signal set rotated 
by an angle of 13.2825° so as to maximize the CPD. For 
STBC-CR the symbols are from a QPSK signal set and rate 1/2 
COD from 16-QAM signal set. For rate 3/4 COD, the symbols 
are chosen from 6-PSK for a throughput of 1.94 bits/sec/Hz 
which is close to 2 bits/sec/Hz. The average transmitted power 
is equal in all the cases i.e. E{tr(S'^S')} = 4, so that average 
energy per bit using the channel model of (|5j is equal. The Fig. 



Observe that rate 3/4 COD and STBC-CR have almost similar 
performance at 2 bits/sec/Hz, and around 1.6 dB coding gain 
over rate 1/2 COD. A possible apparent inconsistency of these 
withtljie results in [32], [33], which report coding gain of over 
is due to the fact that symbol error rate (SER) vs. p 
is plotted in [32], [33]. As rate 1/2 COD chooses symbols 
from 16 QAM and STBC-CR from 4 QAM, SER vs. p plot 
gives an overestimate of the errors for STBC-OD as compared 
to STBC-CR and bit error rate (BER) vs. Et/No is a more 
appropriate plot for comparison at the same through put (2 
bits/sec/Hz). 

From the Table lVll which gives the coding gains of various 
schemes at spectral efficiencies of 2,3,4 bits/sec/Hz, we see 
that the coding gain of STBC-CR and CIOD are nearly equal 
(differ by a factor of 1.11) and significantly greater than 
other schemes. But, the main factor in favor of CIOD as 
compared to STBC-CR (as also any STBC other than STBC- 
OD) is that CIOD allows linear complexity ML decoding 
while STBC-CR has exponential ML decoding complexity. 
At a modest rate of 4 bits/sec/Hz, CIOD requires 64 metric 
computations while STBC-CR requires = 65, 536 metric 
computations. Even the sphere-decoding algorithm is quite 
complex requiring exponential complexity when M < N and 
polynomial otherwise [38]. 

For 4-QAM and 16-QAM constellations, Fig. |5l shows the 
performance for CIOD, STBC-CR and Diagonal Algebraic 
Space Time (DAST) codes of [34]. As expected CIOD shows 
better performance. Finally note that the performance of full- 
diversity QODs [26], [27] is same as the performance of 
CIODs, however QODs are not single-symbol decodable. 
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D. Maximum Mutual Information (MMI) of CIODs 

In this Subsection we analyze the maximum mutual in- 
formation (MMI) that can be attained by GCIOD schemes 
presented in this section. We show that except for the Alamouti 
scheme all other GLCOD have lower MMI than the corre- 
sponding GCIOD. We also compare the MMI of rate-one 
STBC-CR with that of GCIOD to show that GCIOD have 
higher MMI. 

It is very clear from the number of zeros in the transmission 
matrices of GCIODs, presented in the previous sections, that 
these schemes do not achieve capacity. This is because the 
emphasis is on low decoding complexity rather than attain- 
ing capacity. Nevertheless we intend to quantify the loss in 
capacity due to the presence of zeros in GCIODs. 

We first consider the TV = 2, M = 1 CIOD. Equation 0, 
for the CIOD code given in (I40> with power normalization, 
can be written as 



one receiver as is seen from Figures |5] Q and (S) this is 
because the increase in capacity is small from two to four 
transmitters in this case. The capacity loss is substantial when 
the number of receivers is more than one, as these schemes 
achieve capacity that could be attained with half the number 
of transmit antennas. This is because half of the antennas are 
not used during any given frame length. 

Another important aspect is the comparison of MMI of 
CODs for three and four transmit antennas with the capacity 
of CIOD and GCIOD for similar antenna configuration-we 
already know that for two transmit antennas and one receive 
antenna, complex orthogonal designs, (Alamouti code) achieve 
capacity; no code can beat the performance of Alamouti code. 

It is shown in [37] that 



Co(3,Af,p) = ^C(3M,l,^Mp) 



(134) 



V = ^Hs + N 



(130) 



where Co{N,M,p) is the MMI of GLCOD for N transmit 
and M receive antennas at a SNR of p. Similarly, 



where 



^00 









H = 

1^ and where so = so/ 



Co(4,M,p) = ^C(4M,l,^Mp). 



(135) 



and s — [sq sij" , and where so = so/ +JSiq,si = 
sii + jsoQ,so,si G A. If we define Cd{N,M,p) as the 
maximum mutual information of the GCIOD for N transmit 
and M receive antennas at SNR, p, then 

]E{logdetil2 + pH^H)) 



CDi2,l,p) = 



Elog{il + p\hoo\^){l + p\h,o\^)} 
^i^logjl + p\hoo\^} + ii?log{l + p|/iiop} 



C(l,l,p)<C(2,l,p). 



(131) 



It is similarly seen for CIOD code for = 4 given in ( I41> 
that for 



H 



Equation (I135> is plotted for Af = 1, 2 in Fig. Inland ( I133t is 
plotted in Fig. along with the corresponding plots for CIOD 
derived from ( I132> and ( I133t . We see from these plots that the 
capacity of CIOD is just less than the actual capacity when 
there is only one receiver and is considerably greater than 
the capacity of code rate 3/4 complex orthogonal designs for 
four transmitters. When there are two receivers the capacity 
of CIOD is less than the actual capacity but is considerably 
greater than the capacity of code rate 3/4 complex orthogonal 
designs four transmitters. 

Next we present the comparison of GCOD and GCIOD for 
N > A. Consider the MMI of GLCOD of rate K/L. The 
effective channel induced by the GLCOD is given by [37] 
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Cd{^I.p) = 



(^logdet 
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2 ^ 
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where v is a 2K x 1 vector after linear processing of the 
received matrix v, x is a 2K x 1 vector consisting of the 
in-phase and quadrature components of the K indeterminates 
xo, • • • ,xk-i and n is the noise vector with Gaussian iid 
entries with zero mean and variance ||H|p/2. Since ( I136> is a 
scaled AWGN channel with SNR = ||H||2 and rate K/L, 
the average MMI in bits per channel use of GLCOD can be 
written as [37] 



i£;iog{l + |(|/ioop4 

+ lE\og{l+P-{\h2o? 
C(2,l,p)<C(4,l,p) 



f |/»30p 
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Lp I 
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(137) 



)} 



(132) 



observe that H is a iVx A/ matrix. Since ||H|p = H^H where 
H is the NM x 1 vector formed by stacking the columns of 
H, we have 



and 



1 



Co{N,M,p) 



Coii. 1, p) = ■^{C{2, 1, 1, p)} < C(3, 1, p). (133) 

Therefore CIODs do not achieve full channel capacity even 
for one receive antenna. The capacity loss is negligible for 



^CiMN,l,^p) 



(138) 
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where iWi follows from [2, eqn. (10)]. For GCIOD, rec- 
ollect that it consists of two GLCODs, 8i,02 of rate 
K/2Li,K/2L2 as defined in (|79}. Let Ci,o, (^2,0 be the MMI 
of 81, 62 respectively. Then the MMI of GCIOD is given by 

Cd{N,M,p) = j{LiCi^o + L2C2,o} (140) 



{L,Co{N,,M,p) 



+L2Co{N2,M,p)} 



(141) 



2L 

+C I MN2, 1, 



K 
2L2MP 
K 



(142) 



The above result follows from the fact that the GCIOD is block 
diagonal with each block being a GLCOD. When Li = L2 
i.e. 81 = 82 we have 



L~ \ 2 KJ 



(143) 



as we have already seen for = 2, 4. 

Let AC = Cd- Co- For square designs {N ^ L ^ 2°6, b 
odd) we have 



AC = ^C(M2"-i6,l,2"+iMp) 
'-C(M2''6, l,2''Mp). 



N 



(144) 



It is sufficient to consider h = I. When iV = 2, ^ = 2±i = 



1 and AC = C(Af, l,Mp) - C(2M,l,Mp) < 0, as seen 
from [2, Figure 3: and Table 2]. When N > 2, 2a > a + I 
and lima^oo ^ = 2. Also C(A'f2'^-i, l,Mp) is marginally 
smaller than C{M2°',l,Mp) for Af > l,a > 1 as can be 
seen from [2, Figure 3: and Table 2]. It therefore follows that 



Theorem 39: The MMI of square CIOD is greater than 
MMI of square GLCOD except when iV = 2. 

It can be shown that a similar result holds for GCIOD also, 
by carrying out the analysis for each N. We are omitting N ~ 
5, 6, 7. For > 8 we compare rate 2/3 GCIOD with the rate 
1/2 GLCODs. The MMI of rate 1/2 GLCOD is given by 

CoiN,M,p) - ^C{MN,l,2Mp). (145) 
The MMI of rate 2/3 GCIOD is given by, 

CoiN, M, p) = i {C(2Af, 1, Mp) + C{M{N - 2), 1, 2Mp)} 

(146) 

For reasonable values of N that is iV > 8, 

C{MN,l,2Mp) « C{M{N - 2),l,2Afp) and 
C{2M,l,Mp) « C(MN,l,Mp) w C(AfA^, 1, 2A/p) 
and it follows that 



Figure |5] shows the capacity plots for = 8, observe that 
the capacity of rate 2/3 GCIOD is considerably greater than 
that of rate 1/2 GLCOD. At a capacity of 7 bits the gain is 
around 10 dB for M — 8. Similar plots are obtained for all 
N > 8 with increasing coding gains and have been omitted. 
Finally, it is interesting to note that the MMI of QODs is same 
as that of CIODs; however QODs are not SD. 

VII. Single-symbol decodable designs for 

RAPID-FADING CHANNELS 

In this section, we study STBCs for use in rapid-fading 
channels by giving a matrix representation of the multi-antenna 
rapid-fading channels. The emphasis is on finding STBCs 
that are single-symbol decodable in both quasi-static and 
rapid-fading channels, as performance of such STBCs will be 
invariant to channel variations. Unfortunately, we show that 
such a rate 1 design exists for only two transmit antennas. 

We first characterize all linear STBCs that allow single- 
symbol ML decoding when used in rapid-fading channels. 
Then, among these we identify those with full diversity, i.e., 
those with diversity L when the STBC is of size L x N,{L > 
N), where N is the number of transmit antennas and L is 
the time interval. The maximum rate for such a full-diversity, 
single-symbol decodable code is shown to be 2 /L from which 
it follows that rate 1 is possible only for 2 Tx. antennas. The 
co-ordinate interleaved orthogonal design (CIOD) for 2 Tx 
(introduced in section llVt is shown to be one such full-rate, 
full-diversity and single-symbol decodable code. (It turns out 
that Alamouti code is not single-symbol decodable for rapid- 
fading channels.) 

A. Extended Codeword Matrix and the Equivalent Matrix 
Channel 

The inability to write ^ in the matrix form as in (|5jl 
for rapid-fading channels seems to be the reason for scarce 
study of STBCs for use in rapid-fading channels. In this 
section we solve this problem by introducing proper matrix 
representations for the codeword matrix and the channel. In 
what follows we assume that A/ = 1, for simplicity. For a 
rapid-fading channel Q can be written as 



Cd{N, Af, p) « -C(A/iV, 1, 2Mp). 



(147) 



Note that in arriving this approximation we have used the 
property of C(Af , A^, L) that for A^ = 1, as Af increases the 
increment in C is small and also that for a given Af, A", C 
saturates w.rt. p. 



V 



SH + W (148) 

where V € C^^^ (C denotes the complex field) is the received 
signal vector, S <E C^^^^ is the Extended codeword matrix 
(ExCM) (as opposed to codeword matrix S) given by 

' So " 
5i 







Sl-1 



(149) 



where St = [ sq* su ■ ■■ S(^,„i)f J , ff e 
denotes the equivalent channel matrix (EChM) formed by 
stacking the channel vectors for different t i.e. 

Hi 

H = . where Ht 



Hl-1 



hot 
hit 



'■{N~l)t 
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and W e 



->Lxl 



has entries that are Gaussian distributed 



with zero mean and unit variance and also are temporally and 
spatially white. We denote the codeword matrices by boldface 
letters and the ExCMs by normal letters. For example, the 

ExCM S for the Alamouti code, S = 

by 

S 



Xo 



Xl 
Xi} 



IS given 







Xl 









(150) 



Observe that for a linear space-time code, its ExCM S is 
also linear in the indeterminates Xk, k = 0, ■ ■ ■ , K—1 and can 
be written as 5 = Y.k=o^ XkiA2k + XkQA2k+i, where Ak are 
referred to as extended weight matrices to differentiate from 
weight matrices corresponding to the codeword matrix S. 

1) Diversity and Coding gain criteria for rapid-fading 
channels: With the notions of ExCM and EChM developed 
above and the similarity between (|5j and ( I148> we observe 
that, 

1) The distance criterion on the difference of two distinct 
codeword matrices is equivalent to the rank criterion 
for the difference of two distinct ExCM. 

2) The product criterion on the difference of two distinct 
codeword matrices is equivalent to the determinant 
criterion for the difference of two distinct ExCM. 

3) The trace criterion on the difference of two distinct 
codeword matrices derived for quasi-static fading in 
[63] applies to rapid-fading channels also-following the 
observation that tr (S^S) = tr {S"S). 

4) The ML metric Q can again be represented as with 
the code word S replaced by the ExCM, S i.e. 

M{S) = tr{{V - SH)"{V - SH)) . (151) 

This amenability to write the ML decoding metric in 
matrix form for rapid-fading channels (I151> allows the 
results on single-symbol decodable designs of section 
II VI to be applied to rapid-fading channels. 

B. Single-symbol decodable codes 

Substitution of the codeword matrix S by the ExCM, S 
in Theorem ^2 leads to characterization of single-symbol 
decodable STBCs for rapid-fading channels. We have. 

Theorem 40: For a linear STBC in K complex vari- 
ables, whose ExCM is given by, S = J2k=o XkiA2k + 
XkQA2k+i, the ML metric, M(S) defined in ( I151> decom- 
poses as M{S) = Y.kMk{xk) + Mc where Mc = -{K - 
l)tT{V"V), iff 

Af Ai+AfAk = 0,0 <k^l<2K-l. (152) 
Theorem |^ characterizes all linear designs which admit 
single-symbol decoding over rapid-fading channels in terms 
of the extended weight matrices. 

Example 7.1: The Alamouti code is not single-symbol de- 
codable for rapid-fading channels. The extended weight ma- 
trices are 



Ao 

A2 



10 

1 

1 
0-1 



,Ai 



j 

-j 

j 
j 



It is easily checked that the pair Aq , A2 does not satisfy 
equation ( I152> . 

C. Full-diversity, Single-Symbol decodable codes 

In this section we proceed to identify all full-diversity codes 
among single-symbol decodable codes. Recall that for single- 
symbol decodability in quasi-static fading the weight matrices 
have to satisfy (I32> while for rapid-fading the extended weight 
matrices, have to satisfy ( I152> . 

In contrast to quasi-static fading ( I152> is not easily satisfied 
for rapid-fading due to the structure of the equivalent weight 
matrices imposed by the structure of S given in ( I149> . The 
weight matrices Ak are block diagonal of the form ( I149> 
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where A 
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-<lxJV 



In other words even for square 



codeword matrix the equivalent transmission matrix is 
rectangular. For example consider the Alamouti code. 
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is not satisfied as a result we have 
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and hence single-symbol decoding is not possible for the 
Alamouti code over rapid-fading channels. 

The structure of equivalent weight matrices that satisfy 
(I152> is given in Proposition 1411 

Proposition 41: All the matrices Ai that satisfy ( I152> . with 
a specified non-zero matrix Ak in (I153> are of the form 








(0) 











a,.,Af-'^ 



(155) 



where = 0, j Vi. 

Proof The the matrix Ak can satisfy the condition of 

-Af^"A^l\ V<. For a given 



Theorem |40| iff A^!^" A^*^ 



t, A^^^^ A^^ is skew-Hermitian and rank one, it follows that 
^W^f^(t) _ UDU^ where U is unitary and D is diagonal 
with one imaginary entry only. Therefore A^*-* = ijcAl*-* 
where c is a real constant-in fact only the values c = 0, 1 are 
of interest as other values can be normalized to 1, completing 
the proof. ■ 
We give a necessary condition, derived from the rank 
criterion for ExCM, in terms of the extended weight matrices 
Ak for the code to achieve diversity r < L. This necessary 
condition results in ease of characterization. 

Lemma 1: If a linear STBC in K variables, whose ExCM 
is given by, S = Yl,k=o ^ki + XkQA2k+i, achieves 
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diversity r then the matrices A2k^ ^2fc+i together have at least 
r different non-zero rows for every k, Q < k < K — \. 

Proof: This follows from the rank criterion of ExCM 
interpretation of the distance criterion. If, for a given k, 
A-zk, A2k+i together have at less than r different non-zero 
rows then the difference of ExCMs, S — S which differ in Xk 
only, has rank less than r. ■ 
The conditions of Lemma^is only a necessary condition since 
either {xki - Xki) or [xkQ - Xki) may be zero for Xk ^ Xk- 
The sufficient condition is obtained by a slight modification 
of Theorem ^] and is given by 

Corollary 42: A Hnear STBC, S = 
Yl,k=o ^kiA2k+XkQA2k+i where Xk take values from a 
signal set Vk, satisfying the necessary condition of Lemma 
[l] achieves diversity r > iV iff 

1) either A^ Ak is of rank r (r different non-zero rows) 
for all k 

2) or the CPD of A^O. 

Using Lemma n with r = L, we have 

Theorem 43: For rapid-fading channel, the maximum rates 
possible for a full-diversity single-symbol decodable STBC 
using N transmit antennas is 2/L. 

Proof: We have two cases corresponding to the two cases 
of Corollary |^ and we consider them separately. 
Case 1: Ak has L non-zero rows Vk. The number of matrices 
that satisfy Proposition |2] are 2, and the maximal rate is 
R = 1/L. The corresponding STBC is given by its equivalent 
transmission matrix S ~ xqAq, where Aq is of the form given 
in (fT53t . 

Case 2: Ak has less than L non-zero rows for some k. As 
Lemma 1 requires ^2fe,^2fe+i to have L non-zero rows, we 
can assume that A2k has ri non-zero rows and has 
non-overlapping L ~ ri non-zero rows. The number of such 
matrices that satisfy Proposition |^ are 4, and hence the 
maximal rate is i? = 2/L. ■ 
From Theorem |43l it follows that the maximal rate full- 
diversity single-symbol decodable code is given by its ExCM 



S = xoiAo + xoqAi + X11A2 + xiqAz, 
where A2k^ ^2fc+i, A: = 0, 1 are of the form 
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where A, B are of the form given in \\53\ with L = ri and 
L = L — ri respectively. Observe that other STBC's can be 
obtained from the above, by change of variables, multiplication 
by unitary matrices etc. Of interest is the code for L = 2 due 
to its full rate. Setting A = [1 0],i? = [0 1] we have the 
ExCM, 











(158) 



a:i/+ja;oQ 
and the corresponding codeword matrix is 

2:o/+ja;iQ 

a;i/+ja;oQ 

Observe that S is the CIOD of size 2 presented in Section 
II VI Also observe that other full rate STBC's that achieve 



(159) 



,b,ce 



(160) 



full diversity can be achieved from S by performing linear 
operations (not necessarily unitary) on S and/or permutation 
of the real symbols(for each complex symbol there are two real 
symbols). Consequently the most general fuU-diversity single- 
symbol decodable code for = 2 is given by the codeword 
matrix 

Xqi + jxiQ b{xQi + jxiq) 

c{xii+jxQQ) a;i/+ja;oQ 
An immediate consequence is 

Theorem 44: A rate 1 full-diversity single-symbol decod- 
able design for rapid-fading channel exists iff L = = 2. 
Following the results of Section Hvl 

Theorem 45: The CIOD of size 2 is the only STBC that 
achieves full diversity over both quasi-static and rapid-fading 
channels and provides single-symbol decoding. 
Other STBC that achieves full diversity over both quasi- 
static fading channels and provides single-symbol decoding 
are unitarily equivalent to the CIOD for two antennas. Note 
that the CIOD for two antennas dose not have any advantage 
in rapid-fading channels over other SD codes in rapid-fading 
channels. 

Remark 46: Contrast the rates of single-symbol decodable 
codes for quasi-static and rapid-fading channels. From The- 
orem 1^ we have the maximal rate is 2/L for rapid-fading 
channels, while that of square matrix OD [11] is given by 
and that of square FRSDD is given by r'°l^^f/f+' 
respectively. The maximal rate is independent of the number 
of transmit antennas for rapid-fading channels. 

VIII. Discussions 

In this paper we have conducted extensive research on 
STBCs that allow single-symbol decoding in both quasi-static 
and rapid-fading channels. We have characterized all single- 
symbol decodable STBCs, both for quasi-static and rapid- 
fading channels. Further, among the class of single-symbol 
decodable designs, we have characterized a class that can 
achieve full-diversity. 

As a result of this characterization of SD codes for quasi- 
static fading channels, we observe that when there is no 
restriction on the signal set then STBCs from orthogonal 
design (OD) are the only STBCs that are SD and achieve 
full-diversity. But when there is a restriction on the signal set, 
that the co-ordinate product distance is non-zero (CPD 7^ 0), 
then there exists a separate class of codes, which we call 
Full-rank Single-symbol Decodable designs (RFSDD), that 
allows single-symbol decoding and can achieve full-diversity. 
This restriction on the signal set allows for increase in rate 
(symbols/channel use), coding gain and maximum mutual 
information over STBCs from ODs except for two transmit 
antennas. Significantly, rate-one, STBCs from RFSDDs are 
shown to exist for 2, 3, 4 transmit antennas while rate-one 
STBCs ODs exist only for 2 transmit antennas. The maximal 
rates of square RFSDDs were derived and a sub-class of 
RFSDDs called generalized co-ordinate interleaved orthogonal 
designs (GCIOD) were presented and their performance ana- 
lyzed. Construction of fractional rate GCIODs has been dealt 
with thoroughly resulting in construction of various high rate 
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GCIODs. In particular a rate 6/7 GCIOD for iV = 5, 6, rate 
4/5 GCIOD for iV = 7, 8 and rate >2/3 GCIOD for TV > 8 
have been presented. The expansion of signal constellation due 
to co-ordinate interleaving has been brought out. The coding 
gain of GCIOD is linked to a new distance called generalized 
co-ordinate product distance (GCPD) as a consequence the 
coding gain of CIOD is linked to CPD. Both the GCPD and 
the CPD for signal constellations derived from the square 
lattice have been investigated. Simulation results are then 
presented for = 4 to substantiate the theoretical analysis 
and finally the maximum mutual information for GCIOD has 
been derived and compared with GLCOD. It is interesting 
to note that except for N = 2, the GCIOD turns out to be 
superior to GLCOD in terms of rate, coding gain and MMI. 
A significant drawback of GCIOD schemes is that half of the 
antennas are idle, as a result these schemes have higher peak- 
to-average ratio (PAR) compared to the ones using Orthogonal 
Designs. This problem can be solved by pre-multiplying with 
a Hadamard matrix as is done for DAST codes in [34]. This 
pre-multiplication by a Hadamard matrix will not change 
the decoding complexity while more evenly distributing the 
transmitted power across space and time. 

An important contribution of this paper is the novel appli- 
cation of designs to rapid-fading channels, as a result of which 
we find that the CIOD for two transmit antennas is the only 
design that allows single-symbol decoding over both rapid- 
fading and quasi-static channel. It turns out that the single- 
symbol decodability criterion is very restrictive in rapid-fading 
channels and results in constant rate. 

Though we have rigorously pursued single-symbol decod- 
able STBCs and, in particular, square single-symbol decodable 
STBCs, much is left to be desired in non-square STBCs. 
Although non-square STBCs are shown to be useless for rapid- 
fading channels, Su, Xia and Xue-bin-Liang [19], [17] have 
shown for STBCs from ODs in quasi-static channels, that 
higher rates can be obtained from non-square designs. Here 
we list some open problems that were not addressed, or partly 
addressed in this paper. 

• Construction of maximal-rate non-square UFSDDs, RFS- 
DDs. However, the construction of maximal rate non- 
square GLCODs (not GCODs) is itself and open problem 
and any contribution in this direction will greatly enhance 
our understanding of non-square FSDDs. 

• Proof (or contradiction) of existence of non-square FS- 
DDs, 5*, such that 5*^5 is not unitarily-diagonalizable by 
a constant matrix. In Subsection [V] we have shown that 
such square designs do not exist. It would be interesting 
to see if we can obtain even an example of such a design. 
If such a design does not exist then class of UFSDDs 
reduces to GLCODs. In this case the classification of 
UFSDD is complete. Consequently, 

• classification of non-square RFSDD, UFSDD is an open 
problem. In-fact complete classification of RFSDDs ap- 
pears to be even more difficult. Interestingly, [54] shows 
that there exist RFSDDs, that do not belong to the class 
of GCIODs. 

• Even the smaller problem of maximal rates (and design) 
for non-square GCIOD is an open problem. 



• the CPD of non-square lattice constellations and the 
GCPD for both square and non-square lattice constel- 
lations needs to be quantified. It is worth mentioning 
that the authors presented a class of non-square RFSDDs 
called ACIODs in [54] whose coding gain depends on 
CPD and not on GCPD as is the case for GCIODs. 

* Finally, characterization of non-linear STBCs with SD 
property is another open problem. One results in this 
direction is [64]. 

Similar characterization of double-symbol decodable designs 
will be reported in a future paper 

While the final version of the manuscript was under prepa- 
ration the authors became aware of the work [57] that claim 
to unify the results of [48], [54] which is incorrect. The class 
of codes of [57] do not intersect with the class of codes of 
[48] and [54] for the weight matrices of the codes of [57] are 
unitary matrices whereas that of the codes of [48] and [54] 
are not. Furthermore, the STBCs presented in [57], [58] and 
[59] are SD STBCs that do not satisfy (|3T) and such full-rank 
SD STBCs are not considered in this paper. 
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8) Figure |8j The maximum mutual information (ergodic) of CIOD code for four transmitters and one, two receivers 
compared with that of code rate 3/4 complex orthogonal design for four transmitters and the actual channel capacity. 
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TABLE I 

The Encoding And Transmission Sequence For N =2, Rate 1/2 CIOD 





antenna 


antenna 1 


time t 


xoi 





time t + T 








TABLE II 

The Encoding And Transmission Sequence For N =2, Rate 1 CIOD 





antenna 


anteima 1 


time t 







time t + T 





Xll +}X0Q 



TABLE in 

Comparison of rates of known GLCODs and GCIODs for all N 



Tx. Antennas 


GLCODs 


GCIOD (rate-efficient) 


CIOD (delay-efficient) 


N=2 


1 


1 


1 


N=3,4 


3/4 


1 


1 


N=5 


2/3 


6/7 


3/4 


N=6 


2/3 


6/7 


3/4 


N=7 


5/8 


4/5 


3/4 


N=8 


5/8 


4/5 


3/4 


N=2m-3, odd 


(m)/2(m-l) 


2(m+l)/(3m+l) 


7/11 


N=2m-2, even 


(m)/2(m-l) 


2(m+l)/(3m+l) 


7/11 



TABLE rV 

Comparison of delays of known GLCODs and GCIODs N <8 



Tx. Antennas 


GLCODs 


GOOD (rate-efficient) 


GCIOD (delay-efficient) 


N=2 


2 


2 


2 


N=3,4 


4 


4 


4 


N=5 


15 


14 


8 


N=6 


30 


14 


8 


N=7 


56 


35 


8 


N=8 


112 


50 


8 



TABLE V 

The optimal angle of rotation for QPS K and normalized GCPDj^f^ _ for various values of N = Ni + N2. 



N 


Ni 


N2 




Oopt 


GCPDjv,,jv,/4<i^ 


3 


2 


1 


0.555 


29° 


0.3487 




4 


1 


0.5246 


27.76° 


0.28 


5 


3 


2 


0.5751 


29.9° 


0.3869 




4 


2 


0.555 


29.9° 


0.3487 


6 


3 


3 


0.61 


31.7175° 


0.4472 




5 


2 


0.543 


28.51° 


0.3229 


7 


4 


3 


0.5856 


30.35° 


0.40 




7 


2 


0.53 


27.94° 


0.29 


9 


5 


4 


0.591 


30.622° 


0.4135 


10 


8 


2 


0.526 


27.76° 


0.3487 


5 


5 


0.61 


31.7175° 


0.4472 


12 


10 


2 


0.52 


27.5° 


0.265 


6 


6 


0.61 


31.7175° 


0.4472 


N 


N -2 


2 


> 0.5 


> 26.5656° 


> 0.2 
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TABLE VI 

The coding gains of CIOD, STBC-CR, rate 3/4 COD and rate 1/2 COD for 4 tx. antennas and QAM constellations 



R (bits/sec/Hz) 


^CIOD 


^STBC-CR 


^rate 3/4 COD 


^rate 1/2 COD 


2 


0.4478 


0.5 


0.333 


0.2 


3 


0.1491 


0.165 


0.1333 


0.0476 


4 


0.0897 


0.1 




0.0118 
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Fig. 2. Expanded signal sets ^ for ^ = {1, — 1, j, — j} and a rotated version of it. 
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Fig. 3. The plots of CPDi,CPD2 for 9 e [0 90°]. 
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T 



rate 1/2 COD 

rate 3/4 COD 




E,/Nq (dB) 

Fig. 4. The BER performance of coherent QPSK rotated by an angle of 13.2825° (Fig. |2) used by the CIOD scheme for 4 transmit and 1 receive antenna 
compared with STBC-CR, rate 1/2 COD and rate 3/4 COD at a throughout of 2 bits/sec/Hz in Rayleigh fading for the same number of transmit and receive 
antennas. 
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The BER performance of the CIOD with 4- and 16-QAM modulations and comparison with ST-CR and DAST schemes. 
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Fig. 6. The maximum mutual information (ergodic) of CIOD code for two transmitters and one, two receivers compared with that of complex orthogonal 
design (Alamouti scheme) and the actual channel capacity. 
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Fig. 7. The maximum mutual information (ergodic) of GCIOD code for tliree transmitters and one, two receivers compared witli tliat of code rate 3/4 
complex orthogonal design for three transmitters and the actual chaimel capacity. 
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Fig. 8. The maximum mutual information (ergodic) of CIOD code for four transmitters and one, two receivers compared with that of code rate 3/4 complex 
orthogonal design for four transmitters and the actual channel capacity. 
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Fig. 9. The maximum mutual information (average) of rate 2/3 GCIOD code for eiglit transmitters and one, two, four and eiglit receivers compared witli 
that of code rate 1/2 complex orthogonal design for eight transmitters over Rayleigh fading channels. 



